The Burst and Transient Source Experiment (BATSE) 
Earth Occultation Catalog of Low-Energy Gamma-Ray Sources 
Short title: BATSE Earth Occultation Catalog 

B. A. Harmon^ C. A. Wilson, & G. J. Fishman 

NASA Marshall Space Flight Center, SD50, Huntsville, AL 35812 

colleen . a . wilson-hodge@nasa . gov 

V. Connaughton, W. Henze, & W. S. Paciesas 

University of Alabama in Huntsville, Huntsville, AL 35899 

M. H. Finger, M. L. McCoUough? M. Sahi, & B. Peterson 

Universities Space Research Association, SD50, NASA Marshall Space Flight Center, 

Huntsville, AL 35812 

C. R. Shrader 

Universities Space Research Association, Code 661, NASA Goddard Space Flight Center, 

Greenbelt, MD 20111 

J. E. Grindlay 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

D. Barret 

Centre d'Etude Spatiale des Rayonnements, 9 Avenue du Colonel Roche, 
31028 Toulouse Cedex 04, France 

ABSTRACT 

The Burst and Transient Source Experiment (BATSE), aboard the Compton Gamma Ray 
Observatory (CGRO), provided a record of the low-energy gamma-ray sky (^20-1000 keV) be- 
tween 1991 April and 2000 May (9.1y). BATSE monitored the high energy sky using the Earth 
occultation technique (EOT) for point sources whose emission extended for times on the order 
of the CGRO orbital period (^92 min) or greater. Using the EOT to extract flux information, a 
catalog of sources using data from the BATSE large area detectors has been prepared. The first 
part of the catalog consists of results from the all-sky monitoring of 58 sources, mostly Galactic, 
with intrinsic variability on timescales of hours to years. For these sources, we have included ta- 
bles of flux and spectral data, and outburst times for transients. Light curves (or flux histories) 
have been placed on the world wide web. 

We then performed a deep-sampling of these 58 objects, plus a selection of 121 more objects, 
combining data from the entire 9.1y BATSE dataset. Source types considered were primarily 
accreting binaries, but a small number of representative active galaxies, X-ray-emitting stars, 
and supernova remnants were also included. The sample represents a compilation of sources 
monitored and/or discovered with BATSE and other high energy instruments between 1991 and 
2000, known sources taken from the HEAO 1 A-4 (Levine et al. 1984) and Macomb & Gehrels 
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(1999) catalogs. The deep sample results include definite detections of 83 objects and possible 
detections of 36 additional objects. The definite detections spanned three classes of sources: 
accreting black hole and neutron star binaries, active galaxies and supernova remnants. The 
average fluxes measured for the fourth class, the X-ray emitting stars, were below the confidence 
limit for definite detection. 

Flux data for the deep sample are presented in four energy bands: 20-40, 40-70, 70-160, 
and 160-430 keV. The limiting average flux level (9.1 y) for the sample varies from 3.5 to 20 
mCrab {5a) between 20 and 430 keV, depending on systematic error, which in turn is primarily 
dependent on the sky location. To strengthen the credibility of detection of weaker sources (^^5-25 
mCrab), we generated Earth occultation images, searched for periodic behavior using FFT and 
epoch folding methods, and critically evaluated the energy-dependent emission in the four flux 
bands. The deep sample results are intended for guidance in performing future all-sky surveys or 
pointed observations in the hard X-ray and low-energy gamma-ray band, as well as more detailed 
studies with the BATSE EOT. 

Subject headings: gamma rays: observations — methods: data analysis — occupations — surveys 
— X-rays: stars — catalogs 
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1. Introduction 

The Compton Gamma Ray Observatory (CGRO)^, 
one of the four missions in NASA's "Great Observa- 
tories" scries, was launched in 1991 April and oper- 
ated in low Earth orbit until its controlled re-entry 
in 2000 June. CGRO was responsible for many dis- 
coveries in the study of gamma-ray bursts, accreting 
binaries, active galaxies, and pulsars (see general re- 
views by Gehrels & Shrader 1997; Kniffen & Gehrels 
1997; Leonard & Wanjek 2000). The quest for the ori- 
gin of gamma-ray bursts (Meegan ct al. 1992; Fish- 
man & Meegan 1995) led to the development and 
flight of the Burst and Transient Source Experiment 
(BATSE) on the CGRO. In addition to BATSE's pri- 
mary science goals, its nine years of nearly contin- 
uous operation and all-sky capability allowed mon- 
itoring of the low-energy gamma-ray sky using the 
Earth occultation technique (EOT) . Several new ob- 
jects were discovered in the BATSE dataset with the 
EOT, and numerous outbursts of transient and peri- 
odic sources were detected during the mission. Here 
we present a compilation of BATSE EOT observa- 
tions that can be used as the basis for more detailed 
studies on individual sources, for which it is impor- 
tant to have a knowledge of the long-lived emission in 
the BATSE energy range. It is our hope that these 
results will stimulate and guide future observations 
of the low-energy gamma-ray sky. In this work, we 
performed a careful assessment of systematic error in 
deriving results for the entire mission. The intent 
was to push the instrumental and technical limits of 
the present application of the BATSE Earth occulta- 
tion technique (EOT). These data are made available 
on-line to other investigators through the High En- 
ergy Astrophysics Science Archive Research Center 
(HEASARC) at NASA Goddard Space Flight Cen- 
ter, the web site of the BATSE Earth Occultation 
Team at the NASA Marshall Space Flight Center, 
and the Compton Observatory Science Support Cen- 
ter (COSSC). 

To make the presentation of this work tractable. 



^Acronyms and abbreviations used throughout the text and 
tables are listed in Appendix A. 



we divided the catalog into two parts. First, we em- 
phasized the all-sky monitoring capability of BATSE. 
For the bright sources that vary on timescales greater 
than a few hours, we give an overview of activity and 
representative flux and spectral behavior. We also 
briefly make note of the science investigations per- 
formed with the EOT during the mission. Second, 
we performed a "sky survey" of a comprehensive set 
of low-energy gamma-ray sources. We note that this 
is not a true survey, in that it is biased by preselection 
of sources. The emphasis is on obtaining good qual- 
ity results along the Galactic plane where accreting 
sources dominate the low-energy gamma-ray sky. 

We note that other missions with all-sky or wide- 
field coverage were operational between 1991 and 2000, 
A non-exhaustive list includes the Japanese Ginga 
satellite, with its All-Sky Monitor, functioning be- 
tween 1987 February and 1991 October (Tsunemi et 
al. 1989), and covering the bandpass of 1—20 keV, 
the scanning and imaging instruments aboard the 
French-Russian spacecraft Granat (Brandt, Lund, & 
Rao 1990; Roques et al. 1990; Sunyaev et al. 1990; 
Paul et al. 1991), which operated between 1989 De- 
cember and 1997 October, and the U.S. Rossi X- 
Ray Timing Explorer (RXTE) mission, with its All- 
Sky Monitor (2-12 keV) (Bradt, Rothschild, & Swank 
1993; Levine et al. 1996), which began operation in 
1995 December and remains operational. Granat/Sig- 
ma (Roques et al. 1990), a coded mask, sodium io- 
dide imaging detector system, provided considerable 
coverage and detailed images of the central region of 
the Galax;y in the 35—1300 keV bandpass. There 
were also scanning instruments of common design 
on board Granat and Danish EURECA spacecraft 
known as WATCH (Lund 1986; Brandt 1994; Castro- 
Tirado 1994) that monitored a portion of the X-ray 
sky in the 5-120 keV band. The powerful Wide Field 
Cameras ( Jager et al. 1997) on the Italian-Dutch Bep- 
poSax satellite (1.8-28 keV) provided coded aperture 
imaging for both surveying and monitoring of large 
areas of the sky. BATSE complemented these in- 
struments by full-sky monitoring of the energy range 
from 20 keV up to about 1 MeV. Much potential re- 
mains for investigations using combinations of these 
datasets. 



-4- 



BATSE served as an all-sky monitor using the 
Earth occultation technique (Harmon et al. 2002) be- 
tween 1991 April 23 and 2000 May 28 (~9.1y). This 
capability was used to trigger numerous pointed ob- 
servations with the CGRO and other observatories. 
We extracted Earth occultation measurements from 
the BATSE dataset at a maximum rate of twice per 
92m orbit of CGRO in sixteen energy channels be- 
tween 20 kcV and 1 MeV. Intensity data for any part 
of the sky, subject to the geometric constraints of oc- 
cultation, were available for the duration of the mis- 
sion with the exception of a few brief intervals (see 
Appendix B). Although the BATSE EOT was hmited 
in sensitivity by high backgrounds and modest spa- 
tial resolution (-^1°), we could investigate intensity 
and spectral variations over a dynamical range from 
days to years. Although ~daily timescales have been 
accessible with pointed instruments, BATSE was the 
first all-sky monitor since those of Vela 5B (Conner, 
Evans, & BcHan 1969) and Ariel 5 (Holt 1976), to 
cover the upper end of this time range. 

On the other hand, no low-energy gamma-ray 
full-sky surveys have been performed since the pio- 
neering HEAO 1 A-4 experiment in 1977-79 (Levine 

et al. 1984). That experiment covered the energy 
range between 13 and 180 keV, and resulted in the 
detection of about 70 sources (7 of them extragalac- 
tic). Such a survey with BATSE has been attempted 
by our group in the past (Grindlay et al. 1996; Har- 
mon et al. 1997) using the Earth occultation imaging 
technique (Zhang et al. 1993; Harmon et al. 2002). 
This did not prove to be an efficient way to survey 
the sky due to the limitations on computer time com- 
bined with residual features and the small size (a few 
degrees) of the images. A new effort, incorporating 
more sophisticated background modeling techniques 
developed for the International Gamma-Ray Astro- 
physics Laboratory{mTEGRAL){Lei et al. 1999) and 
detector point spread functions for limb geometry, 
holds much promise for achieving this goal (Shaw ct 
al. 2001, 2003). In Fig. 1 we compare the sensitiv- 
ity of the BATSE deep sample (this work) with that 
of the HEAO 1 survey and a proposed future low- 
energy gamma-ray sky survey mission the Energetic 
X-ray Imaging Survey Telescope (EXIST) (Grindlay 



et al. 2001, 2003). The binning of the BATSE sensi- 
tivity curve reflects the maximum energy resolution 
obtainable with the EOT. 

In addition to the list of monitored sources accu- 
mulated during the operation of BATSE, we reviewed 
the HEAO 1 A-4 hard X-ray (Levine et al. 1984) cat- 
alog and the recent gamma-ray source compilation 
of Macomb & Gehrels (1999) for completeness of the 
known high energy Galactic sources, and performed 
a literature search of other high-energy objects above 
the BATSE sensitivity threshold. We did not include 
the EGRET "unidentified" sources from Macomb & 
Gehrels (1999) due to the probable extragalactic na- 
ture of many of them, and the relatively large un- 
certainty in location (~1°) that would strongly affect 
the measured flux using the EOT. 

The sources addressed in this catalog are shown 
on a sky map in Galactic coordinates in Fig. 2. 

We have included 179 sources in our deep sample of 
the sky. The majority of sources are accreting bina- 
ries; however, we include the brighter sources from 
the other object classes for comparison. Of the 179 
sources, most (76%) are accretors, i.e., black hole or 
non-pulsed neutron star binaries (97), or binaries con- 
taining pulsars (41). Again we stress that our deep 
sampling of the known gamma-ray sources should not 
be considered as a true background-limited survey to 
a given brightness threshold. However, by virtue of 
our comparisons to HEAO 1 A-4 survey data, we be- 
lieve that we have a reasonably complete sampling 
of high energy accreting systems in the Milky Way 
Galaxy down to an average flux level of ~30 mCrab 
(20-100 keV). 

We also note that to date, samples of objects 

have been performed with BATSE for active galaxies 
(Connaughton et al. 1999; Malizia et al. 1999; Bassani 
et al. 2000), supernova remnants (McCollough et al. 
1997), globular clusters (Ford et al. 1996), X-ray stars 
(White ct al. 1994), and prompt emission from nearby 
novae (Hernanz et al. 2000). The most promising 
classes of sources, other than accreting binaries, for 
which BATSE has shown positive detections, are ac- 
tive galaxies, such as Seyferts and blazars, and su- 
pernova remnants. As part of our 179-source sam- 
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pie, we included selected supernova remnants (17) 
(some of which contain pulsars or soft gamma-ray 
repeaters (SGRs)), active galaxies (12), and X-ray 
emitting stars (12), that overlap these other studies. 

In Sec. 2, we give a brief overview of BATSE 

and the use of the Earth occultation technique. Wc 
then discuss production details for the catalog. In 
Sec. 3 we present highlights of the all-sky monitor- 
ing effort, and give an overview of science obtained 
for the different classes of detected sources. Sec. 4 
presents the deep sample results and how they are or- 
ganized. We describe a sky-dependent systematic er- 
ror model and apply it to the measured source fluxes. 
We then examine the consistency of source detections 
by focusing on some of the weak sources in the few 
mCrab intensity range. This is done by testing with 
methods that are independent of the occultation step 
fitting routines. Finally, in Sec. 5, we briefly discuss 
applications of the catalog data. 

2. Analysis Methods and Treatment of 
Point Source Data 

2.1. The BATSE Instrument and the Earth 
Occultation Technique 

The Earth occultation analysis technique utilizes 
the large area detectors (LADs) on BATSE, which 
are sensitive to photons above 20 keV. The Compton 
Gamma Ray Observatory including placement of the 
BATSE detector modules on the spacecraft is illus- 
trated pictorially in Fig. 3. The LADs are composed 
of sodium iodide (Nal(Tl)) crystals, 1.27 cm (0.5 in) 
thick by 50.8 cm (20 in) across (2025 cm^ total area 
of one detector). Eight modules are mounted on the 
corners of the CGRO with normal vectors perpen- 
dicular to the faces of a regular octahedron. Any 
point on the sky can viewed by four detectors at an- 
gles less than 90° to the source direction. For the 
EOT, the measurement sensitivity of the LADs is 
maintained by combining statistics from two, three 
and ultimately four detectors at successively larger 
angles from any one of the detector normals. A full 
description of the BATSE detectors and the exper- 



iment can be found, for example, in Fishman et al. 
(1984, 1989a, b). In a companion paper to this catalog 
(Harmon et al. 2002, see also references therein), we 
discussed the use of the EOT in various forms, and 
described in detail a number of aspects of the tech- 
nique such as source identification, sensitivity, and 
systematic error. The Earth occultation technique 
was also developed and applied to BATSE data in a 
separate, parallel effort at the Jet Propulsion Labo- 
ratory (JPL) (Skelton et al. 1994; Ling et al. 1996, 
2000; Ling & Wheaton 2003). We discuss our results 
in comparison to this work in Appendix C. 

In one orbit around the Earth, two occultation 

step features, a rise and set pair, will be superim- 
posed on the background counting rate as each point 
source is occulted. A measurement can be made of 
the intensity of a source in each energy channel at 
rise or set. In practice, two measurements per orbit 
were not always achieved. The most common reasons 
were passages through the lower Van Allen radiation 
belt at the South Atlantic Anomaly when the de- 
tector high voltage was off, or that CGRO was out 
of line-of-sight contact with the NASA Tracking & 
Data Relay Satellites, or data were flagged and not 
available for analysis. In addition, high declination 
sources (| S |>41°) experience an interruption of oc- 
cultations near the orbital poles (see Harmon et al. 
2002), and source confusion, where occultations of 
one source arc indistinguishable from another, limits 
the number of usable occultation steps. The com- 
bined impact of these effects caused Earth occultation 
coverage averaged over one precession cycle (^52d) 
to range between 80%-90%, at best, and at worst, 
about 50%. The count rate for a source in the LAD 
was extracted by simultaneously fitting occultation 
step features with terms for each source in the fit 
and a quadratic polynomial to represent the detector 
background. The fit was performed independently for 
each LAD and in each energy channel. 

Occultation steps (lasting ~10s) are relatively 
sharp features superimposed on the slower background 
variations caused by orbital motion around the Earth. 
We assumed that the background was smooth and 
adequately fit by a second order polynomial. Occa- 
sional failures of this assumption were usually caused 
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by the presence in the data of bright pulsars, weak 
bursts, solar flares, and other disturbances on the 
timescale of the fitting window. Generally, in the 
absence of unusual background fluctuations, we have 
found that the LADs were consistently reproducible 
in the measured fluxes based on tests using the Crab 
Nebula fluxes in a given detector, angle and energy 
bin. The only exception to this was the late mis- 
sion (approx. last two years) loss of gain stabilization 
in LAD 7. We took this considerably further, how- 
ever, to reduce measurement inconsistencies between 
LADs that might be intermittent (unusual background 
fluctuations), and LAD-rcsponse related (time-inde- 
pendent effects). We did this by (1) performing a 
flagging procedure for outliers at the single step level 
(see Harmon ct al. 2002, Sec 2.2, page 152), and also 
outlined in Sec 2.3 of this manuscript), and (2) de- 
veloping an empirical correction to the detector re- 
sponse based on the Crab data in channels 1-10 in 
all eight LADs at angles between and 90 degrees 
(see Harmon et al. 2002, Sec 3.5, pages 167-170). AU 
light curves and spectra in the catalog were subjected 
to these two procedures except for a few sources as 
noted in Sec. 2.3. Any remaining non-Poissonian 
components in the data must then be accounted for 
in estimating systematic errors, which we have done 
on an average basis using a grid-interpolation pro- 
cedure as discussed in Sec 4.2. More details of the 
properties of Earth occultation features and how the 
EOT was used can be found in Harmon et al. (2002). 

2.2. Treatment of Catalog Data 

A raw count rate history for each point source 
was accumulated from fitting occultation steps on 
successive days (usually about 15-30 steps per day). 
The instrument response was deconvolved in a fol- 
lowing step to produce intensities in units of pho- 
ton number, energy flux, or units relative to a refer- 
ence source, such as the Crab Nebula. The process 
of generation and analysis of the occultation count 
rate histories is shown schematically in Fig. 4. The 
raw light curves, uncorrected for detector response, 
were generated with a software package known as 
OF_W. The code carries or collects all information on 



source location, spacecraft orientation and trajectory, 
detector gains, and a limited amount of data on rela- 
tive source intensity, required to extract source count 
rates from occultation step features in the BATSE 
CONT (2.048s, 20-1000 kcV in 16 energy channels) 
data (see Fishman et al. 1989a, b, for a more complete 
description of datatypes). 

In Table 1, we show the standard input param- 
eters for OF_W. The definitions of these parameters 
are: (1) 2t, length of the background data sample 
in seconds. The data that define the fitting window 
were selected to be symmetric about the occultation 
step of the source of interest (SOI). From an anal- 
ysis of relative size of systematic error to the sta- 
tistical error using the second-degree polynomial for 
the background (Shaw et al. 2001), a value for r of 
~110s was found to be optimal. We also allowed an 
additional non-adjustable time of 20s for the width 
of the occultation step itself. (2) 0cut is the maxi- 
mum detector aspect angle of a source (i.e., the angle 
between the detector normal vector and the vector 
toward the source) allowed for source measurement. 
Since there was no improvement in the signal-to-noise 
ratio by combining data from detectors beyond about 
60°, this was a convenient cutoff for choosing which 
detectors were used. This value produced a combina- 
tion of one to four LADs depending on CGRO orien- 
tation. The only exceptions to this were cases where 
a small range of azimuthal viewing angles were oc- 
casionally blocked by the CGRO high gain antenna. 
When this occurred we excluded these LADs from the 
measurement. (3) Oint is the maximum detector as- 
pect angle for a source, with an occultation step time 
in the fitting window of the SOI, that must be simul- 
taneously fit with the SOI (Harmon et al. 2002, see 
for fitting model details) . A value of 70° was selected 
because bright sources such as the Crab or Cygnus X- 
1 still have substantial flux even at very large aspect 
angles. (4) 6 is the minimum time separation between 
occultations of two sources in the fitting window that 
were allowed separate terms in the fit. If the times 
of the SOI and an interfering source fell within the 
fitting window and satisfied the aspect angle criteria, 
but were within (5==10s, then the fit was rejected by 
OF_W. (5) A is the minimum fraction of data packets 
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accepted for an occultation step fit. The acceptance 
criterion helped to minimize potentially singular ma- 
trices and other instabilities in the fitting routines. 
For this work, A was set at 95%. Finally, (6) ^cr is 
the threshold photon flux, which if exceeded, caused 
OF_W to add a fitting term for an occultation step 
of a source within the fitting window of the SOI, pro- 
vided criterion (4) was not invoked. A database was 
kept of sources that were found from previous passes 
through the BATSE dataset to exceed (/>cr for inter- 
vals of more than a few minutes, and typically a day 
or longer. The value of (j)cr, 0.01 photons cm~^s~^, 
refers to the energy band of 20-100 keV, which we 
took as our standard band for this work. 

In addition to a list of sources for measurement, 
OF_W requires the following databases as also noted 
in Fig. 4. 

2.2.1. Flare and Eclipse Databases 

To optimize the sensitivity of the EOT, we tried 
to minimize the number of source terms (SOI -|- po- 
tentially interfering sources) in the fitting window. 
A database of source outburst times, and intensity 
bands provides information to determine whether terms 
for other sources should be included in the fit. This 
database; was built up as new sources were found ei- 
ther through occultation analysis (light curves and 
images) or observations with other instruments. 

Information about source outburst intensity lev- 
els as a function of time was read from the database 
by OF_W before the flux measurement was performed. 
For a first complete pass through the BATSE dataset 
(which we refer to as the "first iteration" ) , (j>cr was set 
at 0.02 photons cm-^s"! (20-100 keV, ~75 mCrab). 
Sources dimmer than this threshold, however, were 
found to introduce significant interference in our mea- 
surements. For this work (the second iteration), we 
lowered (per to 0.01 photons cm~^s~^ (20-100 keV, 
^35 mCrab). In Table 2, we show the sources that 
consistently exceeded 4>cr (identified as "persistent 
sources"), along with their average flux, and for Ta- 
ble 3, we show those that exceeded (per for an iden- 
tifiable, but limited period of time. Although a wide 



variety of source behavior was observed with BATSE 
(See also the assessment of the RXTE ASM data by 
Bradt et al. (2000)), Table 3 sources were broadly 
classified as "transient sources", although duty frac- 
tions could range from <1% up to 90%. If the be- 
ginning and ending times of the fitting window fell 
within the outburst interval of a potentially interfer- 
ing transient source, and it met the usual geometric 
criteria, a fitting term was included. For this partic- 
ular iteration of the catalog, we have added a number 
of low-level outbursts for the transient sources in Ta- 
ble 3 not previously reported in other publications. 
Some special notes regarding the detection and mon- 
itoring of transients are given in Appendix A. 

A few bright X-ray binaries containing neutron 
stars, Her X-1, 4U 1700-377, Vela X-1, Cen X-3, and 
OAO 1657-415, undergo complete eclipses by their 
companion stars and have well-documented ephem- 
erides (see Fig. 5). We took advantage of this in the 
fitting process by eliminating interfering source terms 
for these sources during eclipse. Comparison of the 
eclipse times to the fitting window times were made 
in the solar system barycentric frame (Standish et al. 
1992). 

2.2.2. CGRO Orientation History 

CGRO was 3-axis stabilized, and maintained the 
same orientation for extended periods, typically 2- 
week intervals. The CGRO orientation was defined 
by the right ascension and declination of X and Z- 
axes of the spacecraft in J2000 coordinates, and was 
used to determine the optimum detectors for the oc- 
cultation measurement of the SOI and potentially in- 
terfering sources. 

2.2.3. BATSE CONT Channel Energy Binning 
Lookup Table 

The sixteen channel CONT data had programmable 
energy bins. These were usually kept to a standard 
set of values (see Appendix B), except for occasional 
changes that were made to allow higher resolution 
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measurements of sources with strong spectral cutoffs 
such as pulsars and SGRs. Energy calibration of the 
bin boundaries or "edges" are required to obtain pho- 
ton fluxes. They are also used to compute small elec- 
tronic deadtime corrections in OF_W, which depend 
on the total observed counting rate (source + back- 
ground). 

2.3. Data Flagging and Computation of 
Photon Fluxes 

Before deconvolution of the detector response, 
we flagged (i.e., excluded) any remaining large out- 
liers from the single step raw light curve files. The 
two reasons for doing this were to place a tighter Umit 
on obviously unphysical fitting results, and to exclude 
data contaminated by background fluctuations that 
flagging by daily visual inspection failed to remove 
(Harmon et al. 2002). Both of these situations were 
relatively rare; however, they could produce spuri- 
ous daily flux averages. We settled on a scheme that 
could be implemented without special treatment of 
individual sources, and was based on a filter for vari- 
ations from an average flux level. In this scheme, 
we averaged the count rate from a source for each 
CGRO observing period, and subtracted the aver- 
age from the rates (channels 1-10) over the same pe- 
riod. We then made a cut on the rates at ±3. So- 
to remove outliers. This approach worked well for 
all sources except those that showed strong varia- 
tions of 0.5-1 Crab on timescales of minutes to a few 
hours. The persistent sources that fell into this cate- 
gory were Cyg X-1, Sco X-1, 4U 1700-377, and Vela 
X-1. A few of the brighter, highly variable transient 
sources also exceeded this variability criterion dur- 
ing outburst. Typically, the percentage of rates re- 
jected in this way was less than 5% for the majority of 
sources, and did not significantly affect the measured 
fluxes and uncertainties. However Cyg X-1 could, 
at times, have rejection rates as high as 13%. We 
therefore converted count rate data to photon fluxes 
for these four persistent sources without this flagging 
procedure. For the bright transient sources, we visu- 
ally inspected the raw light curves during the times 
listed in Table 3, and removed flags on apparently 



good data that were set by the automated flagging 
procedure. 

Two methods were developed to obtain flux and 
spectral information from the BATSE occultation count 
rate data. The first method is a fitting procedure 
used commonly in high energy astronomy called "for- 
ward-folding" (Briggs 1995), where we assumed a spec- 
tral model that was a reasonable representation of a 
source flux in the LAD energy range. Then the se- 
lected model was folded through the instrument re- 
sponse (Pendleton et al. 1995) for the illuminated de- 
tectors to determine model count spectra. The count 
data were flt, by varying the model parameters, in a 
chi-square minimization procedure. Finally the best 
fltting model was integrated over the desired energy 
range to determine the photon flux. This method 
avoids a potentially unstable inversion of the BATSE 
response matrices, which have large off-diagonal ele- 
ments due to the LAD's relatively shallow detection 
depth and lack of Compton suppression. It lends it- 
self well to generation of multi-channel spectra, ob- 
taining photon fluxes in a desired energy range, and 
performing broad-band spectral analyses. For the 
spectral analysis results in Tables 2 and 3, we used 
a widely available code known as XSPEC (Arnaud 
1996), and for the flux history generation, a code de- 
veloped locally for specific use with Earth occultation 
data known as HISGEN (see Fig. 4). Both codes use 
forward- folding to estimate photon fluxes. The sec- 
ond method uses the Crab (Harmon et al. 2002, see 
for an example spectrum) as a standard candle by 
obtaining a ratio of the count rate for a given source 
relative to that of the Crab in the same detector ori- 
entation and energy channel. This is convenient for 
computing hardness ratios, epoch folding and inves- 
tigating the deep sample results discussed later. The 
resulting fluxes were less model dependent than those 
using the forward-folding method. 

Four energy bands in Crab-relative units were se- 
lected to report the results of the deep sample in Sec. 
4: 20-40, 40-70, 70-160, and 160-430 keV. The en- 
ergy binning is somewhat constrained by the CONT 
data channel binning scheme (see Appendix B, Table 
Bl). And, although the LADs have some sensitivity 
to fluxes in excess of 430 keV, we chose not to include 
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a higher energy band due to statistics (a relative lack 
of positive detections at higher energy), and the self- 
consistent treatment of systematic error which was 
based on a study only in the 20-100 keV range. Con- 
version factors from the reported four-band Crab rel- 
ative intensities to photon and energy fluxes can be 
found in Appendix B, Table B2. 



3. BATSE All-Sky Monitor Results 

High voltage was first applied to the Large Area 
Detectors on Truncated Julian Date (TJD) 8362, or 
16 April 1991. Final shutdown occurred on TJD 
11691, or 26 May 2000. Data were available between 
TJD 8362 and 11690, except for a few brief peri- 
ods due to orbital reboosts and occasional mission 
anomalies as listed in Appendix B. Light curves and 
representative spectra for the sources considered per- 
sistent or transient by OF_W (a total of 58 objects) 
can be found on the web at 

http: / / gammaray.msfc.nasa.gov/batse/occultation 

in visual as well as tabular form. Table 2 is a list of 
the "persistent (P)" sources known by the EOT flux 
history generator OF_W (see Fig. 4). Likewise, tran- 
sient (T) sources in the BATSE flare database appear 
in Table 3. The tables carry information about source 
type, characteristic spectral data, outburst times, and 
energy flux for these sources in the 20—100 keV band. 

Of the four general classes of sources represented 
in this catalog (x-ray binaries, supernova remnants/pul- 
sars/SGRs, AGNs and x-ray stars), three arc rep- 
resented in Tables 2 and 3. The brightest sources 
in the BATSE database are clearly dominated by 
X-ray binaries (both HMXBs and LMXBs), but we 
also found one supernova remnant / pulsar (the Crab) , 
and three galaxies (Centaurus A, Markarian 501, and 
NGC 4151). No X-ray stars were bright enough to be 
included. In the following sections we will briefly de- 
scribe science investigations of these source classes 
using BATSE EOT results. 



3.1. Black Hole Systems 

Among the persistent black hole candidates, Cyg- 
nus X-1 (e.g., Ling et al. 1997; Paciosas ct al. 1997; 
Zhang et al. 1997a; Brocksopp ct al. 1999b; Zdziarski 
et al. 2002), Cygnus X-3 (McCoUough et al. 1999), 
GX 339-4 (e.g., Harmon et al. 1994; Rubin et al. 1998) 
IE 1740-29 and GRS 1758-258 (e.g., Zhang, Harmon, 
Liang 1997c; Smith et al. 2002) have been studied 
most extensively using the BATSE EOT. Emission 
in the low-energy gamma-ray range is understood to 
be primarily upscattcrcd low-cncrgy photons by the 
inverse Compton process. BATSE studies have il- 
lustrated both spectral and intensity variations as- 
sociated with state transitions (c.g, high (soft), low 
(hard) and off states) of black hole binaries. Such 
state transitions are thought to be driven by changes 
in the mass accretion rate. A prime example of this 
phenomenon is the 1996 hard-soft-hard state transi- 
tion in Cyg X-1, which was monitored in both X-rays 
(2-12 keV) with RXTE/ASM and low-energy gamma- 
rays with BATSE (Zhang et al. 1997a; Esin et al. 
1998; Zdziarski ct al. 2002). These sources have also 
been found to exhibit rather closely correlated ra- 
dio flux variations: IE 1740-29 (Sunyaev et al. 1991; 
Mirabel et al. 1992), GX 339-4 (Fender et al. 1999), 
and Cyg X-3 (McCoUough ct al. 1999). In Fig. 6, we 
show the 9.1y flux histories for four of the more per- 
sistent black hole candidates. It can be seen, that, in 
addition to the periods of low gamma-ray flux (soft 
state) in Cyg X-1 around TJD 9400 and 10270, all of 
these sources show extended periods of low level ac- 
tivity. The low intensity or "quenched" state of Cyg 
X-3, in particular, precedes violent radio flaring and 
jet activity (McCoUough et al. 1999). 

Bright transients (^200 mCrab or more) such 
as GRO J0422-h32 (Shrader et al. 1994; Callanan et 
al. 1995; van der Hooft et al. 1999; Ling & Wheaton 
2003), GRO J1655-40 (Harmon et al. 1995; Zhang et 
al. 1997b), GRS 1915-^105 (Foster et al. 1996; Pa- 
ciesas et al. 1996; Harmon et al. 1997), and XTE 
J1550-564 (Corbel et al. 2001; Wu et al. 2002) have 
been detected at an average rate of about 1 per year. 
In Fig. 7, we show the flux histories of several bright 
black hole transients. Some have conflrmed optical 
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mass functions with lower limits exceeding the theo- 
retical mass of a neutron star. Most have very bright 
high energy spectral tails extending to several hun- 
dred keV, and thus, even without firm mass limits, 
arc considered likely to be black holes. These tran- 
sients reach maximum brightness within a day to a 
few days with highly variable and generally unpre- 
dictable multiple outbursts over an interval of months 
to years following the initial outburst. The overall 
brightness of the secondary outbursts tends toward 
lesser intensity over time (GRS 1915-1-105 being an 
exception). Sometimes there is an apparent return 
to X-ray quiescence between outbursts, e.g., GRO 
J1655-4 (Orosz et al. 1997). Transitions occur be- 
tween high, low, very high and/or intermediate states 
(Corbel et al. 2001), with spectral shapes and lu- 
minosity levels that appear to be the same as the 
persistent black hole systems. Radio spectral den- 
sity measurements and long baseline interferometry 
of black hole transients at times reveal compact ra- 
dio sources, and, with lesser regularity, extended high 
velocity (~ 0.3 — 0.9c) radio-emitting jets (Mirabel 
& Rodriguez 1994; Hjellming & Rupen 1995; Tingay 
et al. 1995). The jet formation is sometimes associ- 
ated with high energy outbursts (Harmon et al. 1995, 
1997). 

3.2. Low Mass X-Ray Binaries with 
Neutron Stars 

Low mass X-ray binaries containing neutron stars 
detected by BATSE fall into two general categories, 
usually known as "atoll" and "Z" sources. The names 
are based on the X-ray spectral behavior (generally 
below the sensitive energy range of the LADs) where 
four energy bands are plotted as paired hardness ra- 
tios on a two-dimensional diagram (a color-color dia- 
gram). As the source intensity varies with time, the 
combined ratios will tend to trace out a Z or atoll- 
like track on the diagram. The physical differences 
between the two NS LMXB categories are not com- 
pletely known. In particular, atoll sources tend to 
emit X-ray bursts, and have spectral states that re- 
semble those of black hole binaries. 



The realization that many LMXB X-ray burst 
sources produce significant emission above a few keV 
came as a result of observations with Granat/Sigma, 
(Barret & Vedrenne 1994) and BATSE (Tavani & 
Barret 1997; Barret 2001). In pulsars, high energy 
radiation is produced by acceleration of particles in 
the strong magnetic fields (>10^^G) near the polar 
regions of the neutron star. The bursters are thought 
to have substantially weaker fields and should be in- 
efficient particle accelerators. In addition, the pres- 
ence of the neutron star surface is thought to gen- 
erate a large flux of soft photons, effectively cool- 
ing surrounding plasma, and suppressing high energy 
emission. Nevertheless, a number of burst sources, 
e.g., KS 1731-260, IE 1724-308 (Barret & Vedrenne 
1994), Aql X-1 (Harmon et al. 1996), 4U 1608-522 
(Zhang et al. 1996b), 4U 0614-i-091 (Ford et al. 1997), 
and GX 354-1-0 (Claret et al. 1994), show significant 
gamma X-ray flux, at times, in excess of 100 keV. 
They are considerably less luminous in the low-energy 
gamma-ray range relative to black hole LMXB sys- 
tems, ~0.05-0.1LBdd (Barret et al. 2000), which is 
probably related to the presence of the neutron star. 

The spectra of LMXB bursters are similar to 

those of black hole systems, typically consisting of 
a soft, disk component with temperatures of a few 
keV and a low-energy gamma-ray tail. The spectra of 
these sources sometimes show a cutoff of 50-60 keV, 
but not always, and may extend to higher energy. 
Comparisons with soft X-ray measurements indicate 
that these sources exhibit anticorrelations between X- 
ray luminosity and spectral hardness (van Paradijs & 
van der Klis 1994; Ford et al. 1997), and thus, as in 
BHCs, spectral behavior of the LMXB NS systems 
should be governed by mass accretion rate. Whether 
or not NS binaries have equivalent spectral states to 
BH binaries is not clear, but these observations sup- 
port the picture of a neutron star surrounded by an 
accretion disk as in BH systems, differing only at 
the innermost regions of the disk. Models involv- 
ing boundary layers between the disk and neutron 
star surface have been devised, although broad-band 
spectral observations do not as yet strictly require 
conditions that are unique to neutron star binaries. 

In Fig. 8, we show examples of persistent and 
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transient LMXBs containing neutron stars detected 
with BATSE. BATSE observations have shown that 
the gamma-ray flux from neutron star LMXBs varies 
over intervals of a few days or longer. Aquila X-1 was 
observed to have extended outbursts of many months 
in the low-energy gamma-ray band (Harmon et al. 
1997) similar to LMXB BHCs. These outbursts were 
correlated with an increase in brightness of its opti- 
cal counterpart V1333 Aquilae. 4U 1608-522 shows 
similar duration outbursts, one of which was bright 
enough to establish that the high energy spectrum 
had a high energy cutoff of about 60 keV (Zhang et 
al. 1996b). GS 1826-238, a source once thought to be 
a black hole candidate, was observed to produce X- 
ray bursts with BeppoSax (Ubertini et al. 1997), and 
therefore is clearly a neutron star system. The source 
reached a level of ~0.01 ph cm~^s~^ around TJD 
10000, and remained roughly a factor of 2 brighter as 
of mission end. 

There is also mounting evidence that the Z-type 

LMXBs, e.g., Sco X-1 (Strickman & Barret 2000; 
D'Amico et al. 2001), GX 17-1-2 (Di Salvo et al. 2000), 
and GX 349-^2 (Di Salvo et al. 2001), exhibit tran- 
sient high energy radiation as do the atoll sources dis- 
cussed above. These spectral tails carry only a few 
percent of the total luminosity of the source (Barret 
2001), and it is unclear whether the high energy emis- 
sion in the Z sources is the same spectral component 
more commonly seen in the atoll sources. So far, the 
presence of such hard tails has not been significantly 
investigated in the BATSE database, although the 
two bright Z-sources Sco X-1 and GX 17+2 were de- 
tected in the deep sample in excess of 160 keV (see 
Sec. 4). 

3.3. Pulsed Sources 

Accreting X-ray pulsars were primarily studied 
with BATSE using pulse frequency and pulsed flux 
measurements (Bildsten et al. 1997; Finger et al. 1999; 
Wilson-Hodge 1999; Wflson et al. 2002, 2003) that in 
general are more sensitive than EOT data. However, 
EOT fluxes been used to obtain pulse fractions, study 
correlations between the spin-up rate and total flux. 



and for the characterization of apastron outbursts. 
EOT light curves for several accreting X-ray pulsars 
are shown in Fig. 9. 

The most common use of EOT measurements 
has been to estimate the pulse fraction, i.e., the ra- 
tio of the pulsed flux to the mean flux (pulsed -|- 
unpulsed) in conjunction with simultaneous or ap- 
proximately simultaneous pulsed flux measurements. 
These measurements, along with orbital parameters, 
can be used to investigate changes in the accretion en- 
vironment, which can depend on a number of system 
properties such as disk-induced torques on the neu- 
tron star and the circumstellar environment of the 
companion. Formal pulsed fractions, using simulta- 
neous data, have been computed for several sources 
including: GRO J1744-28, 25%, 20-40 keV, TJD 10092- 
10098 (Finger et al. 1996a), GX 301-2, - 50%, 20- 
55 keV (Koh et al. 1997), GRS 0834-430, 10-15%, 
20-70 keV, outbursts in 1991 and 1992 (Wilson et 
al. 1997a), EXO 2030-^375, 36(5)%, 30-70 keV, TJD 
9120-9131 (StoUberg et al. 1999), GRO J1008-57, 
67(4)%, 20-70 keV, TJD 9186-9195 (Bildsten et al. 

1997) , A 0535-h262, varied from > 80% at lowest flux 
to 30% at highest flux in giant outburst (Finger, Wil- 
son, & Harmon 1996b; Bildsten et al. 1997), and GRO 
J2058-f42, 36(3)%, 20-70 kcV, TJD 9971-10018 (Wil- 
son et al. 1998). Less rigorously determined pidse 
fractions, with approximately simultaneous data and 
some model-dependency, have been computed for 4U 
1145-619, widely varying from 28(4)-70(20)%, 20-50 
keV, TJD 8678-10742 (Wflson et al. 1994; Wilson- 
Hodge 1999), GS 1843-^00, 7%, 20-50 keV, TJD 10515- 
10516 (Wilson et al. 1997b), and 4U 0115+63, 40- 
60%, 20-50 keV, TJD 11231 (Wilson et al. 1999). 

Simple accretion theory predicts a correlation 
between bolometric flux F and spin-up rate z> in X- 
ray pulsars, described by a power law j> oc F^^"^ . Us- 
ing EOT data, this correlation has been fitted for 
the giant outbursts of GRO J2058+42 (Wilson et al. 

1998) and A0535+262 (Finger, Wilson, & Harmon 
1996b). This correlation indicates the likely presence 
of an accretion disk, because disk accretion is more ef- 
ficient at transferring angular momentum than wind 
accretion. Furthermore, a correlation between the 
frequency of quasi-periodic oscillations and the 20- 
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100 keV flux measured with the EOT was consistent 
with predictions of either the beat frequency or Ke- 
plerian frequency models, confirming the presence of 
an accretion disk in A0535+262 (Finger, Wilson, & 
Harmon 1996b). 

Using dynamic power spectra of BATSE EOT 
and RXTE ASM data, Clarkson et al. (2003) showed 
that the ~ 60-day superorbital period previously known 
in SMC X-1 changed smoothly from 60 days to 45 
days and then returned to its former value on a time- 
scale of approximately 1600 days, varying in a co- 
herent, almost sinusoidal, manner. They interpreted 
this variation as consistent with a radiation-driven 
warping model of the accretion disk. 

Pravdo et al. (1995) and Pravdo & Ghosh (2001) 
folded EOT intensity histories at the 41.498 day or- 
bital period of GX 301-2 from 1993 September - 1994 
July (TJD 9258-9561). They discovered that the 
ncar-pcriastron flare was asymmetric and that a much 
weaker near-apastron flare was also present. In later 
data, after the secular spin-up trend reversed, they 
found that the apastron flare became vanishingly small. 
Laycock, et al. (2003) folded GX 301 2 EOT data 
for the entire BATSE mission and found increasing 
flux levels consistent with orbital phase 0.5, that were 
about 10% of the peak near-periastron fitix. For EXO 
2030-1-375, Laycock, et al. (2003) found evidence for 
an increase in flux near apastron, averaging 30-40% of 
the peak flux in EOT data folded for TJD 8363-9540, 
prior to the change from spin-up to spin-down (Wil- 
son et al. 2002). Direct comparison of folded BATSE 
and RXTE ASM light curves from TJD 10114-10724 
suggested that the apastron emission is softer than 
the periastron emission in EXO 2030-1-375. 

Two X-ray pulsars, discovered with BATSE, were 
initially located using EOT data in conjunction with 
pulsed measurements: GRO J1008-57 (Stollberg et 
al. 1993) and GRO J2058-F42 (Wilson et al. 1995, 
1998). GRO J1849-03 (Zhang et al. 1996a,b), which 
was later found to be identical with GS 1845-024 
(Finger et al. 1999), was located using Earth occul- 
tation imaging (see also Appendix A). 



3.4. Active Galaxies 

Several AGN are sufficiently bright in hard x- 
rays to permit spectral and/or variability studies us- 
ing BATSE Earth occultation data alone. In Fig. 10 
we show intensity histories for four of the brightest 
AGN in the BATSE database. An extensive analy- 
sis of the Seyfert galaxy NGC 4151 by Parsons et al. 
(1998) covered six years of BATSE monitoring data. 
Paciesas et al. (1994) analyzed BATSE observations 
of the radio-loud quasar 3C273 during 1991-1993. Us- 
ing BATSE data, MaHzia et al. (2000) discovered 
hard x-ray emission from the high-redshift quasar 
4C71.07 and presented evidence for long-term source 
variability. Dean et al. (2000) discussed BATSE ob- 
servations of variability in 4C71.07 and two other 
blazars, Mrk 501 and PKS 2005-489. Wheaton et 
al. (1996) reported BATSE spectra of the nearby ra- 
dio galaxy Ccn A from 1991-1992 derived using the 
methodology developed by JPL (Skclton et al. 1994; 
Ling et al. 1996, 2000, see also Appendix C). 

In several cases the BATSE data were useful in 
multi-wavelength studies of AGNs. BATSE observa- 
tions were included in a radio to gamma-ray spectral 
analysis of 3C273 by (Lichti et al. 1995). Courvoisier 
et al. (1999) and Turler et al. (1999) investigated cor- 
relations of BATSE results on 3C273 with optical and 
infrared data during a six-year period. BATSE ob- 
servations of the synchrotron component of the TeV 
blazar Mrk 501 during its 1997 high state were used 
in a multi- wavelength study by Petry et al. (2000). 

Finally, the BATSE light curves were useful as 
context information for observations of AGNs by other 
instnmients. Steinle et al. (1998) used the BATSE 
long-term history of Cen A to characterize the source 
state during COMPTEL observations in 1991-1995. 
Similarly, Haardt et al. (1998) used the BATSE light 
curve in interpreting BeppoS AX ohseivations of 3C273 
in 1996-1997. 



4. BATSE Deep Sample Results 

A comprehensive deep sampling of the low-energy 
gamma-ray sky using the present formulation of the 
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EOT is predicated on the sensitivity curve in Fig. 1. 
The improvement in sensitivity over the HE AO 1 A- 
4 survey is only fully realized if systematic error can 
be minimized. In Harmon et al. (2002) we identi- 
fied systematic errors in the application of the EOT 
to BATSE. Improvements to the detector response 
model, and multiple passes through the BATSE dataset 
with a better picture of the bright sources in the sky 
(Tables 2 and 3) , have allowed us to reduce some con- 
tributers to systematic error. These results encour- 
aged us to combine all the mission data together to 
achieve maximum sky exposure. Later in this section 
we incorporate a systematic error model in assessing 
the detection confidence of the sources in our sam- 
ple. Limited applications of the approach we take 
here can be found in White et al. (1994) and Zhang 
et al. (1996a). 

In Table 4, we present the 179 objects in the 
BATSE deep sample in order of increasing right as- 
cension (a). Firm detections, which include the sys- 
tematic error assessment described later, are indi- 
cated in boldface text. The columns from left to 
right are (1) the BATSE catalog name, (2) an al- 
ternate name recognized by the SIMBAD catalog if 
the BATSE name was not recognized, or other well- 
known name for the source, (3) identification of the 
type of source, if known, given the classes defined in 
Fig. 2. (4) An assigned category (A, B, C, N or 
I, see definitions below) based on the significance of 
detection after, and before (in parentheses) correc- 
tion for systematic error and, (5-6) its celestial coor- 
dinates {a,S) (J2000) and (7-8) Galactic coordinates 
{b,l). The next columns from left to right arc the (9) 
detection significance for the 20-100 keV average pho- 
ton flux averaged over the mission (9.1y) before (in 
parentheses) and after correction for systematic error 
^. and finally, tlu^ flux in mCrab (columns 10-13) of 



■^To compute the 20-100 keV flux for sources for which we 
did not liavc the spectral fitting information given in Tables 
2 and 3, we assumed a power law spectral shape with a pho- 
ton index of -3, typical of accreting systems (an intermediate 
value between spectral slopes seen in black holes and those of 
softer sources such as pulsars). Because the computed 20-100 
keV photon flux is sensitive to the assumed power law spectral 
shape, we only report significance. 



each source in the four energy bands 20-40, 40-70, 70- 
160, and 160-430 keV. The Crab relative fluxes, we 
found, are the most accurate method of displaying 
flux information for the deep sample in the absence 
of precise spectral fitting information for each source. 

It is helpful to classify the sources in our sample 
according to a scheme based on their significance of 
detection and variability, and taking into account the 
systematic eff'ects present in our results. We use the 
following categorization: 

A. (Definite detections) Persistent sources in cat- 
egory A have an average flux (9.1y) that exceeded 
a threshold of 0.01 photons cm^^s in the 20-100 
keV band. Transient sources in category A (denoted 
with a dagger (f) in column (4)) exhibited at least 
one identiflable outburst with an average daily flux 
exceeding a threshold of 0.01 photons cm~^s ~^ in 
the 20-100 keV band. Typical category A sources 
show variability on timescales of at least hours to a 
few days, have a readily identiflable period by EFT, 
or variability that correlates well with behavior in 
other wavelength bands. Because category A source 
identifications are independent of duty fraction, the 
detection significance can vary from a statistical zero 
over 9.1y for brief transients, to very high values for 
bright, persistent sources. 54 sources fall into cate- 
gory A: 21 accreting PSR, 29 BH or NS binaries, 1 
SNR (the Crab Nebula), 3 AGN and no x-ray emit- 
ting stars. 

B. (Definite detections) Sources with a 9.1y-av- 

cragc flux that exceeded lOcr, but with average flux 
<0.01 photons cm-^s in the 20-100 keV band. A 
periodicity search or epoch folding analysis may show 
a period associated with the source. After correction 
for systematic error, 29 sources fall into category B: 6 
accreting PSR, 18 BH or NS binaries, 2 SNR, 3 AGN 
and no x-ray emitting stars. In later sections, we test 
the category B sources in several ways for consistency 
of detection. 

C. (Possible detections) Sources with a 9.1y-av- 
erage flux (20-100 keV) with a signiflcance >3cr and 
<10cr. Variability was not investigated; thus, tran- 
sient activity was not considered and could be present 
in the BATSE dataset. After correction for system- 
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atic error, 36 sources fall into category C: 6 accreting 
PSR, 16 BH or NS binaries, 4 SNR, 1 SGR, 5 AGN 
and 4 x-ray emitting stars. 

N. (Non-detections) Sources with a 9.1y-average 
flux (20-100 keV) of significance between ±3a. Tran- 
sient activity from the sources was not considered, 
and could be present in the BATSE datasct. After 
correction of systematic error, 45 sources fall into cat- 
egory N: 7 accreting PSR, 21 BH or NS binaries, 2 
SNR, 2 SGRs, 4 isolated PSR, 1 AGN and 8 x-ray 
emitting stars. 

I. (Indeterminate) Sources with a 9.1y-avcrage 
flux (20-100 keV) of significance <-3ct, indicating con- 
fusion with nearby sources, or a poorly characterized 
systematic error in the surrounding sky region. The 
category I sources represented, for the most part, our 
inability to separate crowded source regions where 
bright sources known to OF_W were close by. Tran- 
sient activity from category I sources was not con- 
sidered and could be present in the BATSE dataset. 
After correction of systematic error, 15 sources fall 
into category I: 1 accreting PSR, 13 binary BH or 
NS, 1 SGR, and no SNR, isolated PSR, AGN, or x- 
ray emitting stars. Since these sources were indeter- 
minate, we do not report 4-band fluxes for these in 
Table 4. 



4.1. Determination of Systematic Errors 
Due to Sky Location 

In addition to the set of 179 known sources, we 
produced light curves for a grid of 162 test or con- 
trol "source" locations over the entire mission. The 
grid covers the entire Galactic plane, with the highest 
density of test locations in the central bulge region, as 
shown in Fig. 11. The grid has points every 3° along 
the plane between -60° and -|-60° longitude with two 
sets of grid points of the same spacing at -t-6 and - 
6° latitude. Beyond ±60° longitude, there are grid 
points at 6° spacing and 0° latitude. The grid light 
curves were subjected to the same flagging process 
for outliers as performed for the known source set. 
Deconvolution of the detector response for the grid 
points was performed using an assumed power law 



of the spectrum with a photon index of -3. Ideally, 
with no systematic error (a perfect sky model) there 
should be no residual flux. However, a few points 
fall very close (less than a degree) to known bright 
sources, which is desirable, as it allows a more accu- 
rate assessment of systematic error due to the effects 
of flnite spatial resolution. On the other hand, some 
grid points may fall very close to unknown sources 
that are not in our bright source model. 

In Fig. 12, we show the average fluxes for the 
9.1y-long light curves at each grid point as a function 
of Galactic longitude. The samples at -6°, 0°, 6° lat- 
itude are shown together for longitudes between -60° 
and 60°. We also show the standard deviation (in 
sigmas) for the one-day flux averages ^ . Average val- 
ues (absolute) for the fluxes are approximately 0.0015 
photons cm~^s~"'^ (~5 mCrab), with the exception 
of outliers near the Galactic center and points such 
as {l,b) of (-24,0), (42,0). We found no global trend 
in longitude or latitude, but there clearly were re- 
gions where known bright sources are concentrated, 
or unknown point source emission has not been well- 
characterized in our bright source model. Artifacts 
due to subtraction of point source fluxes can be man- 
ifest cither as positive or negative measurements in 
the grid (Harmon et al. 2002), which are allowed by 
the fltting method. In speciflc cases we have found 
that simultaneous fits made to weak sources or empty 
fields with nearby bright sources can sometimes gen- 
erate significantly negative measurements. Unknown 
sources could also produce negative residuals, but 
their actual effect is not separable from the subtrac- 
tion of known sources. 

Wc found that some of the positive grid outliers 
mentioned above fall near sources not known in it- 
eration 2 but revealed very significant flux over at 
9.1y average: 4U 1907-1-097, Cir X-1 and H 1624- 
490. Source confusion in the Galactic center region 
(GCR), as expected, also showed a large residual flux 
from unresolved point source emission. It was not 
possible to attribute the excess GCR flux to spe- 
cific sourc(;s. Wc treat the GCR sources, along with 



^An earlier version of this figure, including the first ~7y of 
data, appears in Harmon et al. (2002). 
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the full source set, using a model for systematic er- 
ror obtained from the grid analysis. However extra 
care should be used in interpreting fluxes from weak 
sources (category B and below) in the few degrees (2- 
4°) surrounding this region. There arc deep images 
of the GCR with Grana^/Sigma (Paul et al. 1991), 
for example, that are better suited for determining 
the average emission characteristics of these sources. 
A similar situation is likely to exist for the Magel- 
lanic Cloud regions, which have not been evaluated 
extensively using BATSE data. 

The second moments, or standard deviations had 
a completely different behavior. Unlike the average 

fluxes, standard deviations were a relatively smooth 
function of sky location. These are plotted in Fig. 
12 in units of Gaussian statistical error for the aver- 
age daily fluxes, i.e., a value of 1 in Fig. 12 means 
there was no contribution from any source of error ex- 
cept Poisson statistics (which, of course, is never the 
case!). The standard deviations were clearly broader 
than predicted for a Gaussian distribution about zero 
flux due to systematic errors {1 a = ~0.01 photons 
cm~^s~^). We found broadening factors of about 
20% (between 70° and 180° longitude and peaking 
near a factor of 2 at about -105°) in excess of normal 
statistics. The standard deviations were largest near 
bright, variable sources such as Vela X-l(/=-96.9°) 
and Cygnus X-1 (/=71.3°). These eff'ects are caused 
by well-known properties of the two sources: the long 
period pulses from Vela X-1, and shot noise from Cyg 
X-1. 



4.2. Application of the Systematic Error 
Model 

Correction factors for the flux average (additive) 
and flux average error (multiplicative) for the 179- 
source set were derived as follows: Using the flux and 
standard deviation for the 162-point grid, we inter- 
polated correction factors both for the flux and its 
associated error corresponding to the locations of the 
known 179 source set from the results shown in Fig 
12. If the known source was within Galactic longi- 
tudes of -60° (300°) and -h60° longitude and within 



±6° latitude, we applied an interpolation based on 
the four nearest neighbors on the grid. Outside of 
the bulge region, -60° longitude and -|-60° longitude, 
we applied a two-point interpolation of the nearest 
neighbors on the grid. For any point sources less than 
±60° longitude, and latitudes between 6° and 10° and 
between -6° and -10°, we also applied a two-point in- 
terpolation. Outside of these regions, well above and 
removed from the Galactic plane, there was no cor- 
rection to the flux average, and a 1.4 multiplicative 
factor was applied to the average error. This value is 
a conservative estimate based on the observed stan- 
dard deviations in Fig. 12. 

The detection significance and relative fluxes in 
Table 4, corrected for systematic error, were then cal- 
culated as 

F - Fcor , . 

where llcor is the corrected significance in column 

9, and F^or and ecor £^re the interpolated correction 
factors derived from the grid analysis, for the 20-100 
keV fiux (F) and error (SF), respectively, and 

Fcor 

Bi,cor = Bi 'p~^i 

where the ^i^cor are the corrected Crab relative 
rates in columns 10-13. The statistical errors for the 
relative rates were also multiplied by the Ccor- 

In Fig. 13, we show the distribution of the 9.1 
yr, 20-100 kcV fiux averages in order of brightness, 
beginning with the brightest persistent source, the 
Crab Nebula. The curve has been computed with and 
without systematic corrections. We also show the as- 
signed category, as determined after the corrections 
were applied. We note that transients, and especially 
those with low duty fractions, may have an average 
flux that is much less than the peak flux, and can 
appear anywhere on the curve. Of the 179 sources, 
about 110-120 sources were above the 3a flux limit. 
This agrees well with the combined number of sources 
(119) in categories A through C. Note that category N 
sources with a slightly negative, non-significant fiux, 
and all of the category I sources, whose fluxes are < 
-3a negative, are not visible in the log plot. The main 
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effect of tlic systematic corrections is to smootfi tlie 
trend of the lower part of the intensity curve, which 
shows a rollover at the low flux end. The physical 
interpretation of the rollover is that we are no longer 
able to distinguish significantly measurable flux due 
to systematic and statistical effects at low source in- 
tensity. We are also prone to missing sources in our 
biased sample such as transients or persistent sources 
with fluxes comparable to those in category C. 

There is a possibility that the rollover of the in- 
tensity curve reflects the true Galactic distribution 
of sources. The depth of our sample in flux is signifi- 
cantly beyond the ~70 sources detected in the HE AO 
1 A-4 survey (Lcvinc et al. 1984), even with BATSE's 
higher energy threshold. For example, HMXBs and 
LMXBs are thought, respectively, to be formed in 
binaries with population I and II progenitor stars, 
and thus have sharply differing Galactic distributions 
(Grimm, Gilfanov, & Sunyaev 2002). However, to 
extrapolate information about source populations in 
our data is difficult due to the missing sources and 
the relatively small size of the sample. A deep, unbi- 
ased survey with an instrument such as EXIST, with 
excellent background rejection and spatial resolution 
characteristics, could produce a much more precisely 
determined extension and shape for the curve shown 
in Fig. 13. 

4.3. Consistency of Deep Sample Detections 

Before concluding this work, it is interesting to 
look at the robustness and consistency of the reported 
detections in Table 4 without going into great detail. 
To do this, we examined sources at fiux levels close 
to the sensitivity limit of the EOT, i.e., just above 
the level where systematic effects dominate. 

We saw from the grid analysis that there were 
positive as well as negative flux offsets at the mea- 
sured control locations, as well as a few significant 
outliers. Ideally, interpolation of the offsets should 
allow us to estimate systematic errors in the point 
source set. However, we have noted some shortcom- 
ings of the approach: there are category I sources 
remaining after the model was applied. Another im- 



perfection is that there are cases where a reasonably 
bright source is omitted from our model (Tables 2 
and 3), and we may erroneously subtract "real" flux 
from the statistically-derived average. 

From the analysis of the grid outliers and in the 
results of Table 4, we found that the HMXB PSR 4U 
1907+097 (Chitnis ct al. 1993), and the LMXBs Cir 
X-1 (laria et al. 2001) and H 1624-490 (Balucinska- 
Church et al. 2000) have a raw significance in the 20- 
100 keV band sufficient to be considered as category 
B, but in failing to identify these sources as bright 
sources, they fall into categories C, C, and N, re- 
spectively, after an overestimated correction. We also 
found reasonable agreement between our uncorrected 
fluxes (i.e., before systematic error corrections were 
made) and the observed low-energy gamma-ray fluxes 
from the literature for these sources. The burster 
IE 1724-307 of the Tcrzan 2 globular cluster is an- 
other example; IE 1724-307 was the only category A 
source found in the deep sample. IE 1724-307 has 
rather strong and variable hard X-ray emission (~30 
mCrab) (Barret & Vedrenne 1994, see also Fig. 8) in 
the 20-100 keV band, and was not considered a bright 
source in the current pass through the data. We did 
not alter the systematic corrections to the flux for the 
four sources in the catalog, but it is likely, pending 
other tests, that these sources should be included in 
the bright source model for a future pass through the 
dataset. 

Next we rigorously examined the category B sub- 
set of sources, a total of 29 objects, which we claimed 
are definite detections. These sources are gathered for 
inspection in Table 5. The (model-dependent) aver- 
age fluxes (9.1y) in the 20-100 keV band arc between 
1.2 and 9.6x10""^ photons cm~^s~^, with detection 
significances, including systematic error, varying be- 
tween 10 and 52a. Unlike Table 4, we list both the 
uncorrected average flux and error, and the interpo- 
lated systematic correction based on the source loca- 
tion with respect to the Galactic plane. 

Eighteen of the 29 sources were within 10° of 
the Galactic plane and thus had corrections to the 
given flux average and error. Correction factors to 
the average flux for the category B sources ranged 
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as high as 150% in the GCR. As indicated in Eq. 
(2), if Fcor > 0, it was subtracted from the statistical 
average. If Fcor < 0, it was added to the statistical 
average. As can be seen in column 4 of Table 4, the 
net effect of the systematic error corrections for the 
deep sample was to render a number of category B 
and C sources as a lower category, and some of the 
category I sources as category N. Of the 29 sources in 
category B, 27 did not require recategorization after 
systematic corrections were applied. 

On the other hand, there were a few cases where 
addition of a large correction factor to the average 
flux of a category C or N (non-detection) source re- 
sulted in a total corrected flux level and significance 
equivalent to a category B source. Two of these ap- 
pear in Table 5: 4U 1820-303, a category B(C) source, 
and 4U 1916-053, a category B(N) source. A pre- 
vious survey for hard X-ray emission in the BATSE 
database (Bloser et al. 1996) reported 2a upper hmits 
of 30 mCrab for both 4U 1820-303 and 4U 1916-053 in 
the 20-100 keV band using 10 days of data. Since that 
time, Bloser et al. (2000a,b) followed up the BATSE 
results with observations of both 4U 1820-303 and 
4U 1916-053 with RXTE, as did Church et al. (1998) 
using BeppoSax. Hard tails were observed in both 
sources. In our deep sample, we report a 20-40 keV 
flux over 9.1y in Table 5 from 4U 1820-303 of 11±0.8 
mCrab (uncorrected) in the 20-40 keV band, and a 
corrected flux of 23±1.1 mCrab. The corrected flux is 
a reasonable factor of 4-5 lower than a recent RXTE 
measurement of the hard X-ray flux by Bloser et al. 
(2000b). 4U 1916-053 showed no signiflcant detec- 
tion (0.9±0.7 mCrab) (uncorrected) in the 20-40 keV 
band, but a corrected flux of 38±1.0 mCrab, reason- 
ably close to the RXTE and BeppoSax results, but 
about a factor of 2 higher. Since it is unlikely that 
the source spends most of the time at flux levels more 
than that observed in the pointed observations, this 
suggests that the magnitude of the systematic cor- 
rection for 4U 1916-053 is over-estimated. The sign 
of 

Fcori i-G., the sense of the correction, though, to 
render a significant positive measurement, is correct. 
Of course, without specific knowledge of variability 
of the hard tails in these objects, it is not possible to 
make a definitive conclusion, but it is possible that 



rather large systematic error corrections could be sus- 
pect in our model. This is not too surprising, since 
we are under sampling a rather complex and rapidly 
changing residual flux as a function of sky location as 
shown in the top frame of Fig 12. We arc also correct- 
ing the 4-band fluxes, which span the energy range 
between 20 and 430 keV, with a correction deduced 
for the 20-100 keV band. 

To investigate the impact of the correction fac- 
tors more generally, we computed the relative change 
in signiflcance for all sources in our sample before 
and after the correction was applied. This is shown 
as a histogram in Fig. 14. The distribution is dom- 
inated by the simple multiplicative correction of the 
second moment (centered at -0.3), but in cases where 
the uncorrected, statistically-derived flux was very 
small relative to the systematic correction to the flux, 
Fcor/ F in Eq. (2) was consequently very large. These 
are represented by the outlying wings of the histogram. 
We therefore made symmetric cuts on the distribu- 
tion of fractional changes in signiflcance about the 
median value at -1.05 on the low end of the distri- 
bution and 0.45 on the high end in Fig. 14. 24 of 
the 179 sources had a fractional change in signifi- 
cance that fell outside of the vertical dashed lines. 
We also found that the absolute value of the signifi- 
cance was increased by application of the systematic 
error correction for 16 of the 24 sources beyond the 
cuts, whereas for most sources, the significance was 
reduced. We therefore decided it was better to quote 
the uncorrected relative Crab rates in the four bands 
for these 24 sources in Table 4, including the two cat- 
egory B sources discussed above. The 4-band fluxes 
with no systematic corrections to the statistically- 
derived fluxes are denoted by the four band rates in 
parenthesis. The judgment here is that the uncor- 
rected 4-band fluxes for these objects are closer to 
the true 9.1y averages. 

In the following sections we illustrate the robust- 
ness of the category B source detections using peri- 
odicity searches and Earth occultation imaging. It 
is important to establish the reliability of the EOT 
method at these low flux levels. For example, detec- 
tion of the black hole candidate LS 5039 (McSwain 
et al. 2001) in the low-energy gamma-ray band has 
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not been reported previously in the literature, and 
4U 1812-12, was reported only for the first time in 
Barret, Olive, & Oosterbroek (2003). 

4.4. Periodicity Search of Light Curves 

X-ray binaries sometimes exhibit periodic or qua- 
si-periodic intensity changes, and in most cases, the 
observed periods correspond to the orbital period of 
the binary (Porb)- Super-orbital periods [Piong) also 
have been observed in accreting neutron star sys- 
tems, and have been attributed to disk precession 
(Wijers & Pringle 1999), and/or periodic reversals 
of the spin direction (Bildsten et al. 1997). Laycock, 
et al. (2003) have investigated periodic behavior in 
HMXBs in the BATSE EOT database. Here we use 
the identification of a period or quasi-period in the 
category B sources to add confidence to the detection 
of the fluxes presented in Table 5. 

Searches for periodic behavior in unevenly sam- 
pled data from the EOT are accomplished using the 
Lomb-Scargle algorithm for Fourier transforms (Scar- 
gle 1982; Home & Baliunas 1986), and epoch fold- 
ing techniques in specific cases. We first performed 
Fourier transforms to create power density spectra 
(PDS) of the 9.1-year one-day average light curves 
to search for periods in the flux variations from all 
sources. 

The PDS were scanned in the frequency range 
between 0.002 d'^ (500 days, 6.6 cycles) and O.Sd'^ 
(the Nyquist frequency of evenly sampled data). A 
simple threshold for reporting of a peak was set at a 
significance, or false alarm probability (FAP) (Press 
et al. 1992), of 0.01 (FAP=1 is roughly the level of 
Poisson noise) for the Lomb periodogram. The FAP 
is related to the PDS peak power z (Press et al. 1992) 
as 

FAP{ z) = 1 - (1 - e-')^ (3) 

or if the probability is <1, FAP{ z) « Me'"". M is 
the number of independent frequencies and is related 
to the total number of flux values N as 

M = Tvf (4) 



The oversampling factor 77 was set at a value of 4, and 
the maximum (high end) frequency k for the PDS was 
set at 2, so that M=AN. It should be noted that the 
FAP is strictly a meaningful quantity for a PDS of 
Gaussian distributed data that is normalized by the 
variance of the dataset. As we have shown earlier for 
EOT light curves, the second assumption is clearly 
violated. Thus the FAPs should only be taken as an 
indication of relative signiflcance. 

In Table 6 we list the search results from the PDS 
analysis of the category A and B sources. We limited 
ourselves to peaks in the PDS that corresponded to 
orbital and super-orbital periods found in the liter- 
ature. The far column indicates the FAP computed 
for the power at the peak frequency. It can be seen 
that the FAP ranges from relative large values ~0.01 
(closer to the Poisson noise level but very significant) 
to very small values for persistent, eclipsing category 
A sources such as Cen X-3, Vela X-1, OAO 1657-415, 
and 4U 1700-377 (see also Fig. 5), as well as some 
of the sources with superperiods such as LMC X-4 
and HER X-1. The reported uncertainty in the cen- 
troid frequency is an approximate error based on the 
oversampling factor and not a formally derived error. 
The quoted error can be reduced by proper fitting of 
the PDS and subtraction of the sinusoidal signal from 
the light curve, and the variance recomputed (Horne 
& Balimias 1986). Since we arc interested only in 
confirming detections of the category B sources, this 
analysis was not performed. 

The only category B sources with detected pe- 
riods using FFTs were SMC X-1, LMC X-4, and 4U 
1538-522. AH are eclipsing HMXBs. Both SMC X-1 
(Wojdowski et al. 1998) and LMC X-4 Laycock, et 
al. (2003) were previously known to show orbital and 
superorbital-related variations in the BATSE data. 
4U 1538-522 was already detected by BATSE via its 
pulsed emission Bildsten ct al. (1997), but this is the 
first reported EOT detection of its orbital period of 
3.728d. 

The FFT analysis was performed only on the 1- 
day resolution light curves with a simple cutoff for 
the FAP and is complicated by low frequency noise, 
CGRO orbital precession, and other effects. For ex- 



- 19 - 



ample, Laycock, et al. (2003) reported a marginal 
result for the orbital period of the HMXB pulsar 4U 
1907+097 (8.3745d), a category C(B) source, that did 
not show up in the PDS search. Since only a few cat- 
egory B sources were detected in the PDS, wc carried 
out epoch folding of the single step (twice per CGRO 
orbit sampling) light curves of at known periods for 
selected category B sources. 

In Fig 15(a), we show 9.1y single step light curves 
of the three category B sources folded at the PorU 
given in column 2 of Table 6 and the light curve of 
4U 1907+097 folded at 8.3745d. In Fig 15(b), we 
show single step folded light curves for sources with 
the reported super-orbital periods in column 3 of Ta- 
ble 6. 4U 1538-522, LMC X-4, and SMC X-1 aU 
exhibit sharp eclipse profiles, as well as prominent 
super-orbital variations for the latter two sources. 
We also epoch-folded single step 4-band rates search- 
ing for reported periods from soft X-ray and opti- 
cal measurements among several of the remaining 
category B sources (some perhaps binary-orbit re- 
lated) in X Per (250d, Delgado-Marti 2001), IE 1145- 
614 (14.335d, Ray & Chakrabarty 2001), GX 349+2 
(21. 7h, Wachter 1997), GX 354+0 (63d, 72d, and 
78d, Kong, Charles, & Kuulkers 1998), 4U 1820-303 
(171d, Chou & Grindlay 2001), LS 5039 (4.117d, Mc- 
Swain et al. 2001), SS 433 (~13d, binary period, 65d, 
jet nodding, ^164d, jet precession, Eikenberry et al. 
2001), 4U 1916-053 (199d, Bloser et al. 2000b), 4U 
2127+119 (17.112h, Ilovaisky et al. 1993) and Cyg X- 
2 (9.84, 40.4, 53.7, 61.3, and 68.8d, Paul, Kitamoto, & 
Makino 2000). Among these 10 sources, no convinc- 
ing periodicities were seen in the epoch folded light 
curves. Given that the average fluxes of the HMXB 
eclipsers 4U 1538-522, LMC X-4, and SMC X-1 are 
comparable to the other category B sources in Ta- 
ble 5, it is somewhat surprising that other binaries, 
which have reported periods in X-rays, did not show 
recognizable periodic variations in the epoch folded 
light curves. 

We point out that several of the category B 
LMXBs are high inclination, compact, ultra-short pe- 
riod binaries known as "dippers". These close bina- 
ries have orbital periods of only minutes to hours. 
The actual dips and/or eclipses are very short, ~few 



minutes, and are observed as absorption-induced changes 
in the soft X-ray flux (-keV) (Church et al. 1998). 
Such short duration absorption dips are not acces- 
sible using the BATSE EOT; however, the reported 
superorbital periods are easily within the timescales 
accessible for epoch folding, yet none were seen. More 
work should be done with these sources, but the ap- 
parent lack of periodicity in the LMXBs (6 out of 6) is 
interesting. Assuming the reported X-ray and opti- 
cal modulations are predominantly associated with 
disk precession as in the HMXBs, such effects do 
not appear to modulate the intensity of the observed 
gamma-ray flux in our category B LMXBs. 

4.5. Earth Occultation Imaging of Category 
B Sources 

Due to the limited success of the periodicity searches 
in the category B light curves, we employed Earth 
occultation imaging to localize the gamma-ray flux 
from these sources. The occultation imaging tech- 
nique (Zhang et al. 1993) allows the construction of 
a 2-D image of a region of the sky to an accuracy 
of 1-2° and has been used to identify new transients 
where the source location was previously unknown 
(e.g., GRO .11655-40, Zhang et al. 1994). For more 
details of the method, see Harmon et al. (2002) and 
references therein. 

Most of the 29 category B sources of the EOT 
catalog are strong enough that we might expect to 
detect them in a 15-day integration ('^6.5cr for a 25 
mCrab source in the 20-100 keV band). A non-detec- 
tion, however, does not imply that the EOT signal 
does not come from this source, because the appli- 
cation and interpretation of the imaging technique is 
complex and depends on more than the strength of 
the source signal and spectral signature. 

Each of the category B sources was examined 
over at least one time interval to look for an ex- 
cess in the image at the source location. A suitable 
limb geometry was chosen by displaying the rising 
and setting limbs over an entire CGRO precession 
period. This period will contain all possible limb ge- 
ometries for a given source and allow the selection 
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of the longest consecutive time interval with near- 
orthogonal rising and setting limbs for the best possi- 
ble spatial resolution. Images were generated for ten 
separate source locations within a 10° by 10° field 
of view. The technique is most efficient for imaging 
sources in the center of the field, but a true signal 
should appear when the source is offset from the cen- 
ter. Offsets up to 1° in each direction were considered 
and images were produced for LAD energy channels 
1 through 3. 

Although most of the category B sources exhib- 
ited at least some variability, the only flux-related 
selection of integration time intervals was made for 
those sources that arc eclipsing binary systems or 
systems with a superorbital period. The RXTE ASM 
light curves (http://xte.mit.edu/ASM_lc.html) allowed 
the selection of time periods when the X-ray emitter 
was not eclipsed by its companion and was not oc- 
culted by a warped accretion disk in those systems 
with a known superorbital period. 

Table 7 summarizes the imaging results. 22 of 
the 29 category B sources yielded at least one cred- 
ible image. The difficulty in obtaining a good limb 
geometry is apparent in some of the sources where the 
contours appear elongated. In others, the noisiness of 
the sky was responsible for the lack of image contours 
associated with the category B source. Figures 16(a- 
d) show some of the noteworthy images for the dif- 
ferent types of astrophysical objects seen by BATSE. 
The AGN 3C 273 has a spectrum hard enough to be 
imaged in the 100-300 keV energy range as seen in 
Fig 16(a). Figures 16(c) and 16(d) are both taken in 
the 50-100 keV energy range and show good source- 
centered images of the supernova remnant PSR 1508- 
59 and the X-ray binary EXO 0748-676 respectively. 
Some source images sufii'er from source confusion. In 
particular 4U 1702-429 and 4U 1705-440 (the latter 
a strong low-energy gamma-ray source in its hard 
state (Barret & Olive 2002)) are close to both each 
other and to the category A source OAO 1657-415 
but separate contours delineate the sources in at least 
some of the grid positions. Another case is that of 
4U 1812-12, which is close to GX 17-1-2, a category 
A source, but has a much harder spectrum. Even 
when it is offset from the center of the image, the 



contours may encompass both sources in channel 1, 
but GX 17-1-2 drops out in channel 2 and the con- 
tours tighten around 4U 1812-12. Although these 
cases show poorly resolved contours, the spatial res- 
olution is sufficient to support the reported detec- 
tions of these sources. The ambiguity often disap- 
pears when the source of interest is in the center of 
the field-of-view, as can be seen in Fig 16(b) where 
the flux from 4U1812-12 (also called Ser X-2) swamps 
that from GX 17-1-2 even in the 20-50 keV range. 



5. Use of the Deep Sample Results 

The BATSE EOT deep sample has interesting 

implications for future studies. The deep sample can 
serve as a guide/comparison tool to an unbiased sur- 
vey with other instruments such as INTEGRAL (Win- 
kler 2001) and Swift (Barthelmy 2000). There also 
remains the possibility of an unbiased survey with 
BATSE (Shaw et al. 2001, 2003). This catalog serves 
as a flrst critical analysis of the entire BASTE dataset, 
and these results are made available on the world 
wide web for potential applications, some perhaps un- 
foreseen. 

The spectral information from BATSE were largely 
unexplored in this work. For example, it was not 
practical for us to present an analysis of the long- 
term spectra in the full 16 channel data that are 
available. However, we made limited comparisons of 
our 4-band results to those of the HEAO 1 survey 
Lcvine et al. (1984), long-term flux averages (^5 yrs) 
of the RXTE ASM compiled by Grimm, Gilfanov, & 
Sunyaev (2002), and to pointed instrument observa- 
tions for AGNs and LMXB neutron stars from the 
BeppoSAX and RXTE. This can be done using the 
conversion factors given in Table A3, and the Crab 
Nebula, which averages about 75 c/s in the 2-12 keV 
ASM band. Allowing for the differing observation 
intervals involved, we find that our 4-band relative 
rates compare well in cross-checks with these results. 
This suggests that the BATSE data could be used at 
a next level to obtain duty fractions for these sources 
in various energy bands. 

We find from the category B analysis that some 
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types of sources such as supernova remnants and AGNs 
show hard photon emission even into the 160-430 keV 
band as shown in Table 5. Accreting pulsars, in con- 
trast, have relatively hard spectra near the BATSE 
threshold and sharp spectral cutoffs at higher energy 
consistent with our understanding of their emission 
mechanisms. Across the ASM (2-12 keV) and BATSE 
(20-40 keV) bands, the fluxes reported for the ac- 
creting pulsars are within a few a agreement for the 
two instruments. The high duty fraction LMXB neu- 
tron stars with hard X-ray emission from pointed- 
instrmncnt observations (Barret et al. 2000) have fluxes 
that arc also comparable with the averages given in 
Table 5. 

We caution users not to over interpret some re- 
sults, such as the sources in category C considered 
as possible detections. Some could indeed be real de- 
tections, but represent the limits of the sensitivity 
of the current version of the BATSE EOT, and the 
fluxes are dominated by the systematic corrections. 
For example, a few X-ray stars (White ct al. 1994, 
see also) appear in this category (none appeared in 
categories A or B and most are in category N), but 
the fluxes are too weak for testing with the meth- 
ods employed in Sec. 4. We stress that the category 
C results definitely can suffer from source confusion 
and other systematic errors that could easily produce 
the low-level observed fluxes. As is, the 4-band rela- 
tive rates in Table 4 for the category C and N sources 
should be used only to derive upper limits on the 9.1y 
averages. 
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Table 1. OF_W Input Parameters 



Parameter 


Description 


Value 
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Length of CONT data fitting window 


220s 
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Aspect angle cutoff for "source-pointed" LADs 


60° 


Oint 


Aspect angle cutoff for interfering sources 


70° 
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Exclusion time for interfering source step 
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A 


Fraction of data packets required in fitting window 


0.95 


V, r 


Photon fiux miiiiimim k-vcl for iiitcrfbriug sourc(>s 


0.01 pli cm -s ' 



Tlux value is for the energy band 20-100 keV with the assumed spectral model for 
the source (see Sec. 4) 



Table 2. Occultation (OF_W) Persistent (P) Sources'* 
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SY2 


20- 




PL 


f) ("2 1 5+0 10"! 


8899-890S 


n 


240 


4 93 /8 


CRAB 


SNR 


20- 

Z(VJ 


-200 


PL 


fi* fo DQ+O 02"! 


1 02fifi-1 0272 

-LU Zj U U J- Vy Zi 1 Z 


X 


732 

. 1 OZ; 


5 80 /4 


\^ X vjr./v J. 


HMXR RHP 


90 




PT, 


fl* /I GQ-UO 04^ 

11 11. . Ot/ZCU I 


1 0000 1 000^ 




. ou / 


^i^i "^9 /7 










OTTB 


1 71 o_|_i 

J. ( J. . O— 1 L Z( . 




2, 




19 11/7 




HMXR RHG? 


90 


900 


PT, 


n /o 01 _Ln 9^^ 

O.U I Z, .0 1 ZCU . ZiU J 


1 0Q40 1 0Q44 


n 


.000 


Zi.OW/ U 


GRS1758-258 


LMXR,RHC 


20- 


-200 


PL 


1.8 (2.72±0.35) 


10838-10842 


0, 


,334 


2.55/5 


GXl+4 


LMXB,PSR 


20- 


-200 


OTTB 


37.2±6.9 


9703-9707 


0. 


,460 


5.74/5 


GX17-I-2 


LMXB,NS 


20- 


■150 


PL 


2.81±0.72 


10007-10012 


0, 


,263 


3.85/4 


GX301-2 


HMXB,PSR 


20- 


-150 


OTTB 


8.0 (7.8±1.1) 


8924-8928 


0, 


,927 


1.06/4 


GX354+0 


LMXB,NS 


20- 


■120 


OTTB 


47.2±26.2 


9888-9892 


0, 


,250 


0.11/3 


HERX-1 


LMXB,PSR 


20- 


■120 


OTTB 


10 (10.4±2.1) 


8965-8969 


0, 


,261 


2.81/3 


LMCX-4 


HMXB,PSR 


20- 


■100 


PL 


3.67±0.75 


9284-9289 


0, 


,187 


3.02/3 


NGC4151 


SYl 


20- 


-150 


PL 


2.4 (1.90±0.30) 


9660-9665 


0. 


,091 


2.82/4 


OA01657-415 


HMXB,PSR 


20- 


-200 


OTTB 


26 (36.5±4.3) 


10600-10603 


0, 


,550 


7.03/5 


SCOX-1 


LMXB,NS 


20- 


■150 


PL 


5.0 (5.38±0.09) 


10001-10006 


1, 


,645 


8.69/4 


SMCX-1 


HMXB,PSR 


20- 


-100 


PL 


2.96±0.66 


11323-11330 


0, 


,065 


2.56/3 


VELAX-1 


HMXB,PSR 


20- 


■100 


OTTB 


15 (11.3±0.8) 


11282-11286 


0, 


,584 


4.00/3 



'^Plots and XSPEC input files for those spectra arc available on http://gammaray.msfc.nasa.gov/batsc/occultation/. 

'"In column 5, parameter values with errors result from spectral fits described in this table. Parameter values without errors 



result from historical spectral fits and were used to generate flux histories. If only a single value is listed, the spectral model for 
flux histories was allowed to default to a power law with photon index = 3.0. For the Crab and Cyg X-1, fl* denotes that the 
power law index was allowed to float when flux histories were generated. 



Table 3. Occultation (OF_W) Transient (T) Sources'' 



Name 


Type 


Outburst Times 


Fit Range 


Model 


Param Value'' 


Time Range of Fit 


Flux (20-100 keV) 


xVdof 






(TJD) 


(keV) 






(TJD) 


(erg cm-^s~^xlO~*) 




Pulsed Transients 


2S1417-624 


PSR 


9600-9708 1 1 54*2-1 1 578 


20-130 


OTTB 


27 n 9 1 +1 31 


9645-9650 


0.400 


7.64/4 


4U0115+634 


PSR 


8363-8375 9483-9523 10040-10072 


20-130 


OTTB 


16.1±3.3 


10042-10046 


0.226 


6.97/4 






10084-10093, 11228-11280 














4U1 145-619 


PSR 


8866-8878 9054-9078 9230-9270 


20-130 


PL 


3.77+0.28 


10356-10361 


0.344 


3.24/4 






9424-9445, 9610-9622. 9805-9811. 


















10186-10198. 10354-10363. 10546-10553 














A0535+262 


PSR 


9044-9057, 9180-9205, 9384-9433, 


20-150 


PL 


3.66±0.31 


9499-9503 


0.241 


6.32/4 






Q4QS-9508 














A1118-616 


PSR 


8616-8678 


20-130 


OTTB 


12.9±2.0 


8619-8624 


0.216 


5.06/4 


EXO2030+375 


PSR 


8385-8395. 8427-8439. 8475-8487, 


20-120 


PL 


1.64±0.71 


8477-8482 


0.092 


3.38/3 






8569-8583 8fi1 9-8fi39 8659-8673 


















8705-8718, 8799-8807. 8842-8858, 


















8890-8906, 8936-8950, 9028-9042, 


















9076-9088. 9114-9137, 9166-9182, 


















9214-9228, 9261-9271. 9487-9498. 


















9580-9593, 9626-9640, 10313-10318, 


















10357-10367, 10404-10412, 10500-10508 


















10998-11017 














GROJ1008-57 


PSR 


9182-9209 


20-100 


OTTB 


25 (22.1±3.0) 


9191-9195 


0.400 


7.85/3 


GROJ1744-28 


PSR 


10052-10203, 10435-10560 


20-150 


PL 


3.62±0.34 


10189-10193 


0.660 


4.90/4 


GRO J 1944+26 


PSR 


11062-11682 


20-120 


OTTB 


12 (5.2±3.1) 


11504-11509 


0.144 


1.72/3 


GROJ2058+42 


PSR 


9977-10016 


20-120 


OTTB 


31.2±5.1 


9991-9996 


0.288 


3.96/3 


GRS0834-430 


PSR 


8375-8440, 8490-8550, 8590-8640, 


20-100 


PL 


3.7 (3.09±0.23) 


8816-8821 


0.264 


6.52/3 






8700-8750, 8810-8850, 8950-8990 


















9100-9140 














GS1843+00 


PSR 


10478-10610 


20-100 


OTTB 


31.2±4.6 


10510-10515 


0.364 


2.70/3 


GS1845-024 


PSR 


8404-8424, 8646-8666, 8886-8908, 


20-150 


PL 


2.53±0.39 


10836-10842 


0.254 


4.19/4 






9130-9150, 9370-9398, 9616-9638, 


















9856-9876, 10098-10122, 10342-10362, 


















10582-10608, 10824-10845, 11067-11084, 


















11307-11328, 11550-11571 














XTEJOlll-7317 


PSR 


11116-11203 


20-100 


OTTB 


12 (13.5±9.1) 


11175-11178 


0.143 


2.76/3 


XTEJ1543-568 


PSR 


11610-11660 


20-100 


OTTB 


15 (16.1±6.1) 


11649-11654 


0.141 


0.86/3 


Nonpulsed Transients 


4U 1543-47 


BHC 


8729-8745, 8782-8843, 8880-9042 


20-200 


PL 


2.7 (2.56±0.11) 


8730-8735 


0.562 


7.28/5 


4U1608-522 


NS 


8415-8602, 8620-8965, 9000-9075, 


20-300 


OTTB 


50 (76.7±15.5) 


11511-11516 


0.179 


7.57/6 



9340-9395, 10165-10280, 10380-10408, 



Table 3 — Continued 



Name 


Type 


Outburst Times 


Fit Range 


Model 


Param Value^ 


Time Range of Fit 


Flux (20-100 keV) 
(erg cm~^s~^ Xl0~*) 


xVdof 






(T.TD) 


(keV) 






(TJD) 








10920-11185, 11300-11600 














4U1630-47 


BHC 


10258-10335, 10842-11455 


20-150 


PL 


2.12±0.44 


10889-10894 


0.132 


0.52/3 


AQLX-1 


NS 


8405-8590, 8720-8783, 9137-9593, 
9910-9950, 10232-10360, 10435-10550, 
11357-11486 


20-200 


PL 


2.48±0.24 


9323-9329 


0.196 


3.81/5 


EXO 1846-031 


BHC 


9302-9313 


20-350 


PL 


2.08±0.11 


9302-9307 


0.243 


4.15/7 


GROJ0422-I-32 


BHC 


8840-9100 


20-500 


U 1 1 ri 


110 (199.9±^8.4) 


8950-8955 


0.356 


12.54/8 


GROJ1655-40 


BHC 


9540-9548, 9560-9580, 9590-9612, 
9644-9712, 9728-9740, 9759-9832, 
9920-9963, 10228-10480, 10510-10555, 
10592-10680 


20-300 


PL 


2.8 (2.55±0.09) 
" 


10370-10375 


0.616 


8.55/6 


GROJ1719-24 


BHC 


9252-9365, 9395-9471, 9502-10106 


20-150 


PL 


2.3 (4.09±0.76) 


9423-9427 


0.173 


2.61/4 


CjKbiUu9-4o 


BHC 


9240-9384 


20-500 


PL 


2.0 (2.42±0.09J 


9271-9276 


0.412 


15.03/8 


GRS1739-278 


BHC 


9950-10600 


20-500 


PL 


1.62±0.40 


10518-10523 


0.097 


6.45/7 


GRS1915+105 


BHC 


8740- (e.o.m.) 


20-350 


PL 


2.8 (2.80±0.10) 


10410-10415 


0.718 


11.84/7 


GS1354-645 


BHC 


10742-10855 


20-300 


OTTB 


88.6±16.0 


10800-10805 


0.239 


9.85/6 


GS1826-238 


NS 


10145- (e.o.m.) 


20-120 


PL 


2.07±0.37 


11499-11505 


0.158 


1.86/3 


GX339-4 


BHC 


8430-8650, 8865-9000, 9336-9506, 
9532-9712, 9753-10820, 11208-11310 


20-400 


OTTB 


60 (117.8±16.4) 


8940-8945 


0.408 


6.99/7 


MKN501 


Blazar 


10549-10552, 10554-10555, 10645-10647, 
10682-10684, 10693-10694 


20-500 


PL 


2.82±0.34 


10549-10552 


0.218 


16.5/8 


NOVA.MUSCAE 


BHC 


8400-8500 


20-600 


PL 


2.5 (1.91±0.05) 


8449-8455 


0.233 


8.89/9 


V4641.SGR 


BHC 


11430-11437 


20-200 


PL 


3.5 (3.23±0.13) 


11430-11437 


0.611 


5.14/5 


XTEJ0421+560 


BHC 


10903-10907 


20-100 


OTTB 


26.9±9.4 


10903-10907 


0.191 


0.23/2 


XTEJ1118+480 


BHC 


11547-11561, 11626- (e.o.m.) 


20-500 


PL 


2.1 (2.09±0.14) 


11551-11556 


0.188 


3.24/8 


XTEJ1550-564 


BHC 


11063-11148, 11200-11271, 11640- (e.o.m.) 


20-500 


OTTB 


101 (108.9±4.9) 


11645-11650 


0.964 


11.61/8 


XTEJ1748-288 


BHC 


10967-10977, 10987-11025 


20-300 


PL 


2.2 (2.88±0.12) 


10969-10974 


0.475 


3.60/6 


XTE J 1859+226 


BHC 


11460-11501 


20-500 


PL 


2.1 (2.80±0.14) 


11470-11475 


0.395 


8.80/8 



^Plots and XSPEC input files for these spectra are available on http://gammaray.msfc.nasa.gov/batse/occultation/. 

^In column 6, parameter values with errors result from spectral fits described in this table. Parameter values without errors result from historical spectral fits 
and were used to generate fiux histories. If only a single value is listed, the spectral model for fiux histories was allowed to default to a power law with photon 
index = 3.0. 



Table 4. BATSE Earth Occultation Deep Sample Results'* 



Name 


Alt Name 


Type 


Category" 


(deg) 


(deg) 


Significance 


20-40 


Flux (r 
40-70 


nCrab)*^'^ 

70-160 


160-430 keV 






HIVIXB PSR 




13.02 


-72.43 


302.86 


-44.69 


4.0 


5.6) 


4 3 + 10 


9 + 08 


5+08 


8 4 + 24 


SMCX-2 


SMC X-2 


HMXB PSR 


Q 


13.64 


-73.68 


302.63 


-43.44 


4. 1 


5.7) 


+ 10 


10 + 8 


6 + 8 


15 + 2 4 


XTEJ0111-T317 


XTE JOlll-732 


H^4XB,PSR 


At 


17.80 


-73.29 


300.97 


-43.76 


1.2( 


1.7) 


2.1 ± 2.0 


3.0 ± 1.7 


4.0 ± 1.8 


0.8 ± 5.2 


SMCX-1 


S^^C X-1 


HMXB, PSR 




19.27 


-73.44 


300.42 


-43.56 


28, 7| 


40. 1) 


24.0 ± 0.8 


9.1 ± 0.7 


-0.6 ± 0.8 


-6.8 ± 2.1 


4U0115+634 




HMXB, PSR 


AT 


19.63 


63.74 


125.93 


1,03 


13. 4( 


18.4) 


9.9 ± 0.9 


4,8 ± 0.7 


1,1 ± 0.7 


3,4 ± 2.1 


TTARI 


TT ARI 


Nova-like star 


N 


31.72 


15.30 


148.53 


-43.80 


-1.2 


-1.7) 


-0.4 ± 0.7 


-0.6 ± 0.7 


-2.2 ± 0.7 


-4.2 ± 1.8 


LSI+61_303 




HMXB, Radio Star 


C 


40.13 


61.23 


135.68 


1.09 


5.2 


9.5) 


3.0 ± 0.9 


2.3 ± 0.7 


3.3 ± 0.7 


5.1 ± 2.1 


NGC1275 




SY 




49.95 


41.51 


150.58 


-13.26 


10. 7( 


15.0) 


7.4 ± 0.8 


3.0 ± 0.7 


3.4 ± 0.7 


4.0 ± 2.0 


V0332-|-530 


V 0332-|-53 


t± Jvi A rS , Jr D rt 




53.75 


53.17 


146.05 


-2.20 


1.0 


2.3) 


1.1 ± 0.7 


-0.1 ± 0.7 


-0.1 ± 0.7 


1.5 ± 1.9 


HR.1099 




RS CVn 


r- 


54.20 


0.59 


184.91 


-41.57 


4.6 


6.4) 


2.3 ± 0.7 


1.7 ± 0.7 


0.4 ± 0.7 


2.0 ± 1.8 


XPER 


X PER 


HMXB, PSR 




58.85 


31.05 


163.08 


-17.14 


23. 8( 


33.3) 


9.6 ± 0.8 


12.9 ± 0.7 


6.6 ± 0.7 


0.0 ± 2.0 


VW^HYI 


VW HYI 


DN 


N 


62.30 


-71.29 


284.89 


-38.13 


2.0 


2.8) 


1.3 ± 1.0 


1.0 ± 0.8 


-1.1 ± 0.8 


-8.0 ± 2.2 


XTEJ0421+560 


CI Cam 


HMXBjBHC? 


AT 


64.90 


56.00 


149.17 


4,12 


1.8( 


4.5) 


1.6 ± 0.8 


0.3 ± 0.7 


-0.1 ± 0,7 


1,5 ± 1.9 


GROJ0422+32 




LMXB,BHC 


At 


65.43 


32.91 


165.88 


-11.91 


120. 7( 


169.0) 


37.5 ± 0.8 


45.2 ± 0.7 


56,8 ± 0.7 


57,2 ± 2.0 


HD32918 




RS CVn? 


C 


74.58 


-75.28 


287.41 


-33.24 


3.8 


5.4) 


2.5 ± 1.0 


2.3 ± 0.8 


1.1 ± 0.8 


-5.4 ± 2.4 


SAXJ0501 + 11 


SAX J0501. 7+1146 


LMXB?,NS 


C 


75.43 


11.78 


188.91 


-17.95 


4.4 


6.1) 


2.5 ± 0.8 


1.5 ± 0.7 


1.8 ± 0.7 


-0.5 ± 2.0 


4U0512-041 


4U 0513-40 


LMXB,NS 


N(C) 


78.53 


-40.04 


244.51 


-35.04 


2.9 


4.1) 


0.7 ± 0.8 


2.4 ± 0.7 


0.9 ± 0.8 


2.3 ± 2.1 


SGR0526-66 




SGR 


N 


81.50 


-66.08 


276.09 


-33.25 


1.2 


1.6) 


0.8 ± 1.4 


0.7 ± 1.3 


-0.5 ± 1.3 


-7.6 ± 3.5 


LMCX-4 


LMC X-4 


HMXB, PSR 


B 


83.21 


-66.37 


276.34 


-32.53 


25. 9( 


36.3) 


20.4 ± 1.0 


9.6 ± 0.8 


0.8 ± 0.8 


-8.0 ± 2.2 


CRAB 




S1MR,PSR 


A 


83.63 


22.01 


184.56 


-5.78 


1560.6 


2055.7) 


1021.7 ± 0.8 


1022.5 ± 0.7 


1022.5 ± 0.7 


1038.1 ± 1.9 


SN1987A 




SNR. 


C 


83,87 


-69.27 


279,71 


-31 .94 


4.3 


6.1) 


3.9 ± 1.1 


2.2 ± 1.0 


1.8 ± 1.0 


-6.6 ± 2.7 


A n f^, Q f^, Q 

J\ U o vj- u u ft 








83,92 


-66.86 


276.88 


-32, 20 


-1.5 


-2.0) 


-0.3 ± 2.0 


-l..> ± 1,8 


-2.9 ± 1,8 


-10.2 ± 5.0 


A0535+262 


lA 0535+262 


HMXB,PSR 


At 


84.73 


26.32 


181.45 


-2.64 


46. 6( 


64.7) 


45.1 ± 1.0 


23.6 i 0.8 


7.7 ± 0.8 


^ ^"5 


MCG8-1 1-1 1 


MCG +08-11-011 


SYl 


C(B) 


88.73 


46.44 


165.73 


10.41 


8.6( 


12.1) 


4.5 ± 0.8 


4.6 ± 0.8 


5.4 ± 0.8 


2.5 ± 2.2| 


4U0614-}-091 




LK/IXB^N'S 


B 


94.28 


9.14 


200.88 


-3.36 


40. 4( 


54.0) 


22.6 ± 0.9 


18.1 ± 0.7 


13.5 ± 0.7 


8.2 ± 2.^ 


A0620-00 


lA 0620-00 


LMXBjBHC 


C 














( 0.1 ± 0.7) 


( 3.0 ± 0.6) 


( 1.1 ± 0.6) 


( 7.3 ± li§5 


ijlljJVillN KjA 






C(N) 


98 48 


17 77 


195 14 


4 27 


6 6 


2.7) 


( 0.2 i 0.6) 


( 1.4 + 0.6) 


( 1.8 + 0.6) 


( 3.1 ± 1.5) 


PSR0656+14 




PSR 


N 


104.95 


14.24 


201. 11 


8.26 


1.6{ 


-2.2) 


( -1.0 ± 0.6) 


( -1.0 ± 0.5) 


( 1.0 ± 0.5) 


( 0.7 ± 1.4) 


EXO0748-676 




LMXB,NS 


B 


117,14 


-67,75 


279.98 


-19,81 


15, 4l 


21.6) 


9.6 ± 1.0 


9.1 ± 0.8 


4.7 ± 0.8 


-5.7 ± 2.2 


CPPUP 


CP PUP 


Nova 


N(I) 


122,94 


-35.35 


252,93 


-0.84 


-1.2 




-1.4 ± 1,2 


-0.5 ± 1,0 


-0,1 ± 1,0 


-0.3 ± 3,1 


VELA 




SNR, PSR 


B 


128.84 


-45,18 


263.55 


-2.79 


14. 7( 


24.3) 


10.3 ± 1.2 


10.1 ± 1.1 


11.6 ± 1.1 


12.4 ± 2.8 


GRS0S34-430 


Ginga 0834-43 


HMXB, PSR 


At 


129,16 


-43,26 


262. 16 


-1,45 


25, 8| 


46. 1) 


22.1 ± 1.1 


14.8 ± 1.1 


7.1 ± 1.1 


4.3 ± 3.0 


4U0836-429 


4U 0836-42 


LMXB.NS 


C(B) 


129.35 


-42.90 


261 .96 


-1.12 


7.6( 


12.4) 


5.8 ± 1.3 


6.1 ± 1.1 


6.1 ± 1.1 


3.3 ± 3.2 


OJ287 


OJ 287 


BL Lac 


C 


133.70 


20. 11 


206.81 


35.82 


4.3 


6.1) 


3.2 ± 0.7 


1.1 ± 0.7 


0.7 ± 0.7 


0.6 ± 1.8 


VELAX-1 


VELA X-1 


HMXB, PSR 


A 


135.53 


-40.55 


263.06 


3.93 


250. 6< 


449.0) 


227.2 ± 1.1 


84.1 ± 0.9 


8.1 ± 0.9 


1.7 ± 2.7 


H0918-549 




LMXB,NS 


N(C) 


140. 11 


-55.21 


275.85 


-3.84 


2.3 


9.2) 


1.7 ± 1.2 


1.5 ± 1.0 


0.8 ± 1.0 


-2.5 ± 2.8 


SN1993J 




SNR 


C 


148.85 


69.02 


142. 15 


40.91 


4.8 


6.8) 


2.0 ± 1.0 


2.9 ± 0.8 


3.2 ± 1.0 


8.5 ± 2.4 


GROJ1008-57 




HMXB, PSR 


At 


152.50 


-58,31 


283 .03 


-1.82 


0.3( 


5,1) 


0.3 ± 1.2 


0,2 ± 1.0 


0,0 ± 1.0 


-0.5 ± 2,7 


GRS 1009-45 


Granat 1009-45 


LMXB BHC 


At 


153.40 


-45.08 


275.88 


9.34 


9.9( 


24.1) 










1E1024-57 


2E 1024.0-5733 


HMXB? 




156.49 


-57.81 


284.52 


-0.24 


0.7 


5.0) 


6 ±12 


4 ± 10 


±10 


^0^9 ± 2 7 


1E1048-593 


2E 1048.1-5937 


HMXB? 


N(C) 






288.26 








0.0 ± 1.1 


0.0 ± 1.0 


0.0 ± 1.1 


0.0 ± 2.8 


o rv X u o o- o 






N(C) 


164.49 


-52.45 


285.99 


6.65 


0.5 


^4.1)^ 










MKN421 




BL Lac 


(J 


166.11 


38.21 


179.83 


65.03 


4.4 


6.2) 


2 6 ± 8 


17 ± 7 


2 ± 7 


2 8 ± 2 


XTEJ1118+480 




LMXB, BHC 


At 


169.57 


48.05 


157.62 


62.32 


6.5( 


9.1) 


3.0 ± 0.8 


3.2 ± 0.8 


3.8 ± 0.8 


6.1 ± 2.2 


AlllS-616 


lA 1118-616 


HMXB, PSR 


At 


170.24 


-61.92 


292.50 


-0.89 


-1.6 


2.1) 


( 2.3 ± 0.9) 


( -1.1 ± 0.8) 


( -4.6 ± 0.9) 


( -5.5 ± 2.3) 


CENX-3 


CEN X-3 


HMXB, PSR 


A 


170.31 


-60.62 


292.09 


0.34 


33. 9( 


62.3) 


32.2 ± 1.1 


7.1 ± 1.0 


-2.7 ± 1.0 


-6.8 ± 2,6 


NOVA.MUSCAE 


Nova Mus 1991 


LMXB, BHC 


At 


171.61 


-68.68 


295.30 


-7.07 


6.2( 


13.1) 


4.1 ± 1.1 


4.1 ± 1.0 


3.1 ± 1.1 


0.2 ± 2.9 


AX1138-68 




HMXB, PSR 


C(N) 


174.51 


-69.90 


296.66 


-7.93 


4.8 


0.6) 


( 0.9 ± 0.7) 


( -0.2 ± 0.6) 


( -0.7 ± 0.7) 


( -5.1 ± 1.8) 


1E1145-614 


IE 1145.0-6140 


HMXB, PSR 


B 


176.87 


-61.95 


295.49 


-0.01 


22. 7( 


37.3) 


20.3 ± 1.3 


14.6 ± 1.1 


4.1 ± 1.1 


-5.3 ± 3.1 


4U1145-619 




HMXB, PSR 


At 


177.00 


-62,21 


295.61 


-0,24 


24. 8( 


39.8) 


22.2 ± 1.1 


16,4 ± 1.1 


4,7 ± 1.1 


-5.2 ± 2.9 


SN1994W 




SNR 


N 


180.55 


62.14 


132.70 


54.04 


1.9 


2.7) 


1.6 ± 1.0 


0.8 ± 0.8 


0.9 ± 0.8 


8.5 ± 2.4 


NGC4151 




SYl 


A 


182.64 


39.40 


155.08 


75.06 


93. 0( 


130.2) 


38.6 ± 0.8 


42.2 ± 0.7 


37,0 ± 0.7 


20.9 ± 2.0 


GX301-2 




HMXB, PSR 


A 


186.66 


-62.77 


300.10 


-0,03 


155. 8( 


221.3) 


147.2 ± 1.1 


37.2 ± 0.9 


0,1 ± 0.9 


-3.3 ± 2.6 


3C273 




QSR 


B 


187.28 


2.05 


289.95 


64.36 


51.71 


72.4) 


18.6 ± 0.8 


19.3 ± 0.7 


23,2 ± 0.7 


32.6 ± 2.0 


REJ1255+266 




uvs 


C 


193.78 


26.69 


5.00 


89.07 


4.5 


6.3) 


2.8 ± 0.8 


2.0 ± 0.7 


1.0 ± 0.8 


0.1 ± 2.1 


3C279 




QSR 


C 


194.05 


-5.79 


305.11 


57.06 


5.0 


6.9) 


2.8 ± 0.8 


2.1 ± 0.7 


3.0 ± 0.7 


7.0 ± 2.1 


4U1254-690 


4U 1254-69 


LMXB,NS 


N(C) 


194.40 


-69.29 


303.48 


-6.42 


2.1 


8.4) 


2.0 ± 1.1 


1.2 ± 0.9 


0.0 ± 1.1 


-3.1 ± 2.8 


GX304-1 




LMXB,NS 


I 


195.32 


-61.60 


304.10 


1.25 


-3.7 


-3.8) 










PSR1259-63 




HMXB?, PSR 


N 


195.70 


-63.84 


304.19 


-0.99 


-2.5 


-1.7) 


-3.2 ± 1.1 


-1.2 ± 1.0 


0.0 ± 1.1 


-13.4 ± 2.9 


CEN_A 




SY2 


A 


201.37 


-43.02 


309.52 


19.42 


110. 7( 


154.9) 


40.7 ± 0.8 


46.0 ± 0.8 


54.5 ± 0.8 


66.3 ± 2.2 
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Name 


Alt Name 


Type 


Category^ 




(a 

(d 


eg) 


(dcg) 


Significance 


20-40 


Flux (n 

40-70 


nCrab)*^'^ 

70-160 


160-430 kcV 


4U1323-619 


4U 1323-62 


L]VIXB,1MS 


B 


201, 


,65 


-62.14 


307.03 


0.46 


12. 0( 27.5) 


8.2 ± 1.2 


8.4 ± 1.0 


4.7 ± 1.0 


-0.1 ± 2.8 


GS1354-645 


Ginga 1354-64 


LMXBjBHC 


At 


209. 


,50 


-64.70 


309.97 


-2.74 


4.0( 14.4) 


2.5 ± 1.2 


2.5 ± 1.0 


2.4 ± 1.0 


-0.7 ± 2.8 


NGC5548 




SYl 


B 


214, 


,50 


25.13 


31.95 


70.50 


13. 3( 18.6) 


6.0 ± 0.8 


6.3 ± 0.7 


5.7 ± 0.7 


6.4 ± 2.0 


2S1417-624 




HMXB,PSR 


At 


215. 


,30 


-62.70 


313.02 


-1.60 


15. 5( 26.2) 


12.6 ± 1.2 


9.2 ± 1.0 


1.0 ± 1.0 


-4.5 ± 2.8 


CENX-4 


CEN X-4 


LMXB.NS 


N(C) 


224. 


.59 


-31.67 


332.24 


23.88 


2-2( 3-1) 


0.9 ± 0.8 


0.9 ± 0.7 


1.2 it 0.7 


0.6 ± 2.1 


PSR,1509-58 




SNR,PSR 


B 


228. 


,48 


-59.14 


320.32 


-1.16 


11. 8( 26.5) 


6.9 ± 1.1 


6.9 ± 0.9 


8.2 it 0.9 


9.7 ± 2.4 


H1517+656 




BL Lac 


N 


229. 


.45 


65.42 


102.26 


45.39 


1.5( 2.1) 


0.2 ± 1.0 


1.2 ± 0.8 


1.9 it 1.0 


4.8 ± 2.5 


CIRX-1 


GIR X-1 


LMXB.NS 


C(B) 


230. 


.17 


-57,17 


322.12 


0,04 


3.8( 16-1) 


3,5 it 1.1 


1.2 ± 0.9 


1.2 it 0.9 


0.3 ± 2.6 


TRAX-1 


TRA X-1 


LMXB,NS 


N 


232. 


.07 


-61. SS 


320.32 


-4.43 


-1.1( -1,4) 


-3.4 ± 1.1 


0.2 ± 0.9 


1.2 it 0.9 


-1.1 ± 2.6 


4U1538-522 


4U 1538-52 


HMXBjPSR 


B 


235. 


,60 


-52.39 


327.42 


2.16 


18. 3( 24.6) 


16.3 ± 1.1 


7.7 ± 0.9 


1.3 ± 0.9 


-4.1 ± 2.6 


XTE J 1543-568 




HMXBjPSR 


At 


236. 


,05 


-56.74 


324.99 


-1.46 


-3.7( -6.9) 


-3.0 ± 1.1 


-2.1 ± 0.9 


-3.4 ± 0.9 


-9.1 ± 2.6 


4U1543-47 




LMXB,BHC 


At 


236. 


,79 


-47.67 


330.92 


5.43 


10. 1( 16.6) 


6.4 ± 1.0 


5.0 ± 0.9 


4.3 ± 0.9 


1.5 ± 2.5 


XTEJ1550-564 




LMXB,BHC 


At 


237, 


,74 


-56.48 


325.88 


-1 .83 


17. 1( 23.1) 


10.5 ± 1.0 


12.2 ± 0.9 


7.1 ± 1.0 


-2.3 ± 2.5 


2S1553-342 




HMXB.PyR 


I 


239. 


.46 


-54,41 


327.95 


-0,86 


-5.3( -16,9) 










4U1608-522 




L.IvIXB,NS 


At 


243. 


,17 


-52.43 






50. 9( 43.8) 


21.7 ± 1.0 


31.9 ± 0.9 


29.3 ± 1.0 


9.9 ± 2.5 


SCOX-1 


SCO X-1 


LMXB.NS 


A 


244. 


,98 


-15.64 




23 78 


460. 6( 644.9) 


364.3 ± 0.8 


45.0 ± 0.8 


12.2 ± 0.8 


9.9 ± 2.2 


H1624-490 




LMXB,NS 


N(B) 


247, 


.02 


-49.20 


334.92 


-0.27 


0.0( 24.7) 


0.0 ± 1.2 


0.0 ± 1.0 


0.0 ± 1.2 


0.0 ± 3.0 


4U1627-673 


4U 1626-67 


LMXB,PSR 


B 


248. 


,07 


-67.46 


321.79 


-13.09 


14. 8( 20.7) 


14.5 ± 1.0 


4.2 ± 1.0 


-1.2 ± 1.0 


-7.2 ± 2.7 


4U1630-47 




LMXB,BHC 


At 


248. 


,50 


-47.39 






-10. 2( 48.0) 


( 31.4 ± 0.8) 


( 24.3 ± 0.7) 


( 13.1 ± 0.8) 


( 10.0 ± 2.1) 


SGR1627-41 




SGR 


1(B) 


248, 


.97 


-47.60 


336 97 


^^2 


-5.4( 19.5) 










4U1636-536 




LMXB,NS 


B 


250. 


,23 


-53.75 


332.92 


4 82 


10. 9( 21.3) 


11.9 ± 1.3 


5.2 ± 1.2 


0.6 ± 1.2 


-3.4 ± 3.0 


GX340+0 




LMXB,NS 


c 


251, 


.45 


-45.61 


339.59 


■ 


4.8( 4.8) 


10,8 ± 1.3 


-4.8 ± 1.2 


-0.9 ± 1.2 


-5.8 ± 3.3 


]V1KN501 






At 


253. 


,47 


39.76 


63.60 


38.86 


8.1( 11.3) 


5.4 ± 1.0 


4.4 ± 0.8 


4.9 ± 0.8 


2.7 ± 2.5 


/^TD*^T"i«KK Art 




l_j IVl XJ , ±J XX 


At 








344.98 


2.46 


Ar\ r\\ 


/Ifl A -L "1 A 
<l\J . O i X . O 


35.0 ^ 1.3 


26.9 i 1.5 








J_j iVl 1_> , 1 o XV 




5 


Aft 

■ 


35.3 


58.15 


37.52 




41.7 ± 0.8 


14.1 ± 0.7 


2.5 ±0.8 


1.4 ± 2.1 


O AO 16 57-41 5 




HIVIXB , PSR 








' oo 


344.35 


0.31 


55. 5 ( 62 .6) 














ij i\ 1 yv JlJ J IN O 


N(I) 






, .. 


353.83 


7.27 


-0.4( -12,3) 


- . ± .3 


-0. ± . 


- . ± . 


- . ± . . 


GX339-4 




LMXBjBHC 


A 1 


255. 


,71 


-48.79 


338.94 


-4.33 


89, 0( 116.8) 


58.8 ± 1.4 


71.8 ± 1.2 


73.7 ± 1.3 


58.0 ± 3.5 


4U 1700-377 




H1VIXB,NS? 


A 


255. 


,99 


-37.84 


347.76 


2.17 


193. 0( 260.4) 


212.8 ± 1.6 


149.6 ± 1.3 


63.6 ± 1.5 


13.4 ± 3.9 


GX349-I-2 




LMXBjNS 


B 


256. 


,43 


-36.42 


349.10 


2.75 


14. 4( 17.1) 


21.9 ± 1.6 


5.3 ± 1.5 


2.3 ± 1.6 


0.6 ± 4.2 


401702-428 




LiKlXB ,NS 




256. 


,56 


-43.04 


343.89 


-1.32 


12. 3( 10.7) 


20.2 ± 2.5 


17.0 ± 2.2 


14.1 ± 2.3 


13.4 ± 6.4 


H1705-250 




r TlffVD TT 

1j JVl Js. ±5 , ±5 rllj 


C(I) 








358.59 


9.06 




( -4.0 i 1.0) 


( -3.7 i 0.8) 


( -2.9 ± 0.9) 


( 2.7 i 2.3) 


4U 1705-440 


4U 1705-44 


LflV/IXB NS 




257 


23 


44 10 


343.32 


-2.34 




22.8 i 1.6 






16.5 i 3.9 


PSR1706-44 




PSR 


N 


257. 


.43 


-44.49 


343. 10 


-2.69 


2. 1 ( -0.0) 


( -1.1 zb 3.3) 


( -5.2 zb 3.0) 


( 3.7 it 3.3) 


( 10.2 it 9.1) 


MXB 1715-321 


IH 1715-321 


LMXB NB 


I 


259. 


70 


-32. 18 


354-13 


3.07 


-6.8( - 12.4) 










G R,0 J1719-24 


■ 


LIVIXB ,BHC 


At 


259, 


,90 


-24.97 


0.18 


7.02 


70. 4 ( 78.5) 


38.1 _^ 1.5 


4- 

58.7 ± 1.3 


57.6 ± 1.4 


45 .3 it 3.5 


IE 1724-307 


4U 1724-307 


LIVIXB ,NS 










356.32 


2.30 
















Ij IVl A. iz> , IN G 




202 


93 


16 96 


8.51 


9.04 


12 2( 11 2) 


15 3 ^ 15 


^5^± 1^3* 


^3 4 ± 13 


^7 9 ± 3 e 


GX354-|-0 




LIvlXBjNS 




263. 


,00 


-33.84 


354.30 


-0.16 


50. 2( 68.2) 


46.4 i 1.3 


34.5 ± 1.1 


18.2 ± 1.1 


2.0 i 3.2 


GXl-|-4 




H]\4XB PSR 










1.94 


4.80 


85 ■8( 109.9) 


82.9 i 1.6 








1730-335 


4U 1730-333 


LMXB NS 


N 


263, 


.35 


-33.39 






-1.0( -0.7) 


( -6.7 ± 3.1) 


( 4.2 ± 2.7) 


( 4.4 ± 2.9) 


( 2.8 ± 7.7) 


KS1731-260 




LMXBjNS 


N(C) 


263, 


.55 


-26.09 


1.07 


3.66 


1.6( 6.1) 


4.1 ± 3.3 


2.8 ± 2.8 


0.3 ± 3.0 


-0.1 ± 8.1 


SLX1732-304 




LMXB NS 


I(N) 


263, 


.95 


-30.48 


357.56 


0.99 


-6.3( -0.5) 










GX359+2 




LMXB,NS 


1(C) 


264. 


.65 


-28.48 


359.57 


1.55 


-12. 0( 9.1) 










4U1735-444 




LMXB,NS 


C 


264. 


.74 


-44.45 


346.06 


-6.99 


4.0( 3,1) 


10.3 ± 1.5 


-3.9 ± 1.3 


1.3 ± 1.4 


9.0 ± 3.7 


GRS1737-31 


Granat 1737-31 


LMXB,BHC 


I 


264. 


.95 


-30.96 


357.62 


0.01 


-12. 1[ -4.0) 










XTEJ1739-285 




XRS 


1(C) 


264. 


.98 


-28.48 


359.73 


1.31 


-16. 9( 6.0) 










GRS 1739-278 


Granat 1739-278 


LMXB,BHC 


At 


265. 


,67 


-27.75 


0.67 


1.18 


0.3( 20.3) 


0.6 ± 3.7 


0.5 ± 3.1 


0.3 ± 3.2 


0.2 ± 8.7 


1 El 740-29 


1E1740. 7-2942 


LMXB?,BHC 


A 


265. 


,97 


-29.74 


359.12 


-0.10 


92. 4( 183.9) 


76.7 ± 1.5 


74.3 ± 1.4 


74.6 ± 1.4 


61.2 ± 3.7 


GR.OJ1744-28 




LIVIXB, PSR 


At 


266, 


,14 


-28.74 


0.05 


0.30 


5.5( 27.7) 


17.2 ± 4.3 


11.4 ± 3.8 


4.3 it 4.1 


4.8 ± 10.9 


KS17.'i 1-293 




LMXB,NS 


N{A) 


266. 


.20 


-29.35 


359.56 


-0.07 


0.8[ 14.7) 


3.2 ± 6-8 


2.9 it 6-0 


1.0 it 6.4 


1.3 ± 17.5 


H1743-322 




LMXB,NS 


N(C) 


266. 


.26 


-32.23 


357.13 


-1.61 


2.5( 3,9) 


1.6 ± 2-6 


4,6 ± 2-3 


5.9 it 2-3 


9.4 ± 6-2 


SLX1744-3nO 




LMXB,NB 


C(A) 


266. 


.86 


-30.04 


359.26 


-0.91 


7.7( 44.4) 


10.6 ± 1,9 


7.5 it 1-7 


4.9 it 1,S 


5.2 ± 4-7 


GX3 + 1 




LMXB,1SJS 


B 


266, 


,98 


-26.56 


2.30 


0.79 


13. 1( 29.0) 


15.2 ± 1.4 


8.1 ± 1.3 


5.5 m 1.3 


8.5 ± 3.6 


XTEJ1748-288 




LIvIXB,BHC 


At 


267. 


,03 


-28.48 


0.68 


-0.23 


-12. 4( 49.9) 


( 40.5 ± 1.1) 


( 28.9 ± 1.0) 


( 17.8 ± 1.0) 


( 18.9 ± 2.6) 


H1745-203 




LMXB.NS 


I 


267. 


.22 


-20.37 


7.73 


3.80 


-7.3( -10.3) 










XB01745-248 


GPS 1745-248 


LMXB,NS 


C(B) 


267, 


,23 


-24.89 


3.84 


1.46 


10. 0( 10.8) 


11.6 ± 1.5 


5.4 ± 1.3 


5.2 ± 1.3 


7.5 ± 3.5 


GROJ1750-27 




HMXB.PSR 


N(B) 


267, 


.50 


-27.00 


2.16 


0.17 


1.2( 20.1) 


1.4 ± 1.4 


0.6 ± 1.3 


0.5 ± 1.3 


1.0 ± 3.3 


4U1746-370 


4U 1746-37 


LMXB,NS 


N 


267, 


.55 


-37.05 


353.53 


-5.00 


0.5( 2.7) 


0.9 ± 1.5 


0.5 ± 1.3 


0.2 ± 1.3 


0.3 ± 3.7 


EX01747-214 




LMXB,NS 


I 


267, 


.61 


-21.42 


7.00 


2.95 


-10. 7( -11.4) 











CO 
CO 
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Name 


Alt Name 


Type 


Category^ 




(a 
(d 


eg) 


(dcg) 


Significance 


20-40 


Flux ( 
40-70 


TiCrab)*^'^ 

70-160 


160-430 keV 


SAXJ1750-29 


SAX J1750. 8-2900 


LMXB?,NS 


1(B) 


267 


70 


-29.00 


0.53 


-1.00 


-12. 4( 30.9) 










XTEJl 755-324 




LMXB,BHC 


C(B) 


268 


87 


-32.48 


358.04 


-3.63 


7.1{ 15.0) 


7.7 ± 1.5 


5.5 ± 1.2 


3.8 ± 1.4 


5.9 ± 3.4 


H1755-338 




LMXB ,NS 


C 


269 


67 


-33.81 






8.2( 8.1) 


8.3 ± 1.4 


6.1 lb 1.2 


4.3 ± 1.4 


7.2 ± 3.4 






T "N/rVD "NTC? 










^5^08^ 


1 02 


3.8( 6.0) 


21.2 i 6.8 






-24.0 i 17.5 






1_j1V1 .^±3 f l3 XI 




270 


30 


25 74 


4.51 


-1.36 






42 4 ^^1^1 


50 4 ± 1^ 


OB.U It iJ.U 


GX9+1 




LMXB.NS 


N 


270 


38 


-20.53 


9.08 


1.15 


2.1( 0.2) 


( 6.3 ± 1.0) 


( -5.5 ± 0.9) 


( -7.3 ± 0.9) 


( -8.5 ± 2.4) 


2S1803-245 




LMXB,BHC 


I 


271 


71 


-24.58 


6.15 


-1.90 


-5.9( -3.1) 










SGR1806-20 


AX 1805.7-2025 


SGR 


N 


272 


17 


-20.41 


10.00 


-0.24 


-1.6( -0.5) 


( 1.1 ± 0.9) 


( -1.2 ± 0.8) 


( -2.9 zb 0.8) 


( -9.5 zb 2.1) 


NOVA_SAG 


Nova Sgr 1991 


Nova 


N(I) 


273 


56 


-32.21 


0.20 


-6.97 


2.2( -13-1) 


( -7.5 ± 0.9) 


( -7.2 ± 0.8) 


( -5.5 zb 0.9) 


( -6.2 zb 2.2) 


GX13+1 




LMXB,PSR 


C 


273 


63 


-17.16 


13.52 


0.11 


7.7( 4.0) 


( 1.6 ± 1.0) 


( 3.1 ± 0.9) 


( 0.8 ± 1.0) 


( -2.2 zb 2.5) 


4U1812-12 




L1VIXB,NS 


B 


273 


81 


-12.03 


18.11 


2.40 


16. 2( 22.9) 


22.5 ± 2.5 


25.2 ± 2.2 


21.1 ± 2.3 


15.6 zb 6.0 


GX17+2 




LMXBjNS? 


A 


274 


01 


-14.04 


16.43 


1.28 


65. 4( 85.2) 


56.1 ± 1.2 


32.1 ± 1.0 


24.4 ± 1.0 


12.0 zb 2.8 


AM_HER 




CV 


N 


274 


05 


49.87 


77.87 


25.88 


1.5( 2.1) 


2.3 ± 1.0 


0.4 ± 0.8 


-1.2 ± 1.0 


-2.3 ± 2.4 


V4641_SGR 




LMXBjBHC 


At 


274 


84 


-25.41 


6.78 


-4.79 


-6.9( -4.7) 


( -0.6 ± 0.9) 


( -3.3 ± 0.7) 


( -4.6 ± 0.8) 


( -8.8 ± 2.0) 


4U 1820-303 




LMXB,NS 


B(C) 


275 


92 


-30.36 


2.79 


-7.91 


13. 9( 9.1) 


( 10.9 ± 0.8) 


( -1.5 ± 0.7) 


( -3.5 ± 0.7) 


( -7.9 ± 1.9) 


H1822-000 




LMXB.NS 


N(I) 


276 


34 


-0.01 


29.94 


5,79 


-0,8( -9.9) 


-0.6 ± 1.2 


-0.4 ± 1.0 


-0.5 zb 1.0 


-0.1 zb 2,6 


H1822-371 




LMXB, IMS 


B 


276 


45 


-37.10 


356.85 


-11.29 


17. 7( 24.8) 


21.5 ± 1.1 


2.5 ± 1.0 


-1.9 zb 1.0 


-6.3 ± 2.7 


LS5039 




HMXB,BHC 


B(C) 


276 


56 


-14.85 


16.88 


-1.29 


10. 7( 7.3) 


( 3.4 ± 0.8) 


( 3.0 ± 0.7) 


( 2.0 zb 0.7) 


( 3.0 zb 1.8) 


GS1826-238 


Ginga 1826-238 


LMXB, MS 


At 


277 


37 


-23.80 


9.27 


-6.09 


25. 5( 64.0) 


15.7 ± 1.1 


17.0 ± 1.0 


12.7 zb 1.0 


5.3 zb 2.7 


SCTX-1 


SCT X-1 


HMXB.PSR 


N(I) 


279 


18 


-7. 58 


24.51 


-0.20 


-1.1( -3.3) 


-1.1 ± 1.2 


-0.6 ± 1.0 


0.6 zb 1.0 


2.0 zb 2.6 


SERX-1 


SER X-1 


LMXB,NS 


C(B) 


279 


99 


5.04 


36.12 


4.84 


6.0( 14.2) 


5.0 ± 1.1 


2.7 ± 0.9 


2.2 ± 1.1 


2.5 ± 2.7 


GS1843+00 


Ginga 1843+00 


HMXB,PSR 


At 


281 


50 


1.00 


33.20 


1.67 


18. 8( 36.1) 


11.6 ± 1.0 


11.0 ± 1.0 


7.6 ± 1.0 


6.8 ± 2.6 


GS1845-024 


2S 1845-024 


HMXB,PSR 


At 


282 


07 


-2.42 


30.42 


-0.40 


2.3( 35.9) 


1.4 ± 1.0 


1.4 ± 0.9 


0.8 ± 1.0 


0.5 zb 2.6 


EXO1846-031 




LMXB,BHC 


At 


282 


30 


-3.10 






-1.7( 24.5) 


( 10.3 ± 0.8) 


( 10.9 ± 0.7) 


( 6.7 zb 0.7) 


( 4.6 zb 1.8) 


4U1850-087 


4U 1850-08 


LMXB.NS 


C(B) 


283 


27 


-8.71 


25.36 


-4.32 


8.8( 16.3) 


5.6 ± 1.0 


5.7 ± 0.9 


4.1 ± 1.0 


1.3 zb 2.5 




!XlTEJ1859-j-226 




LMXB BHC 


At 












6.2( -3.6) 


( -2.2 i 0.7) 


( -0.7 i 0.6) 




{ -"K 4.-1-1 R"! 
V 31 J.. 07 ^ 


SGR1900+14 




SGR 


C(B) 


286 


81 


9.32 


43 02 


077 


9.2( 52.1) 


7.2 ± 1.1 


4.8 ± 0.9 


3.0 ± 1.1 


1.1 ± 2.6 


4U1905+000 


4U 1857+01 


LMXB,NS 


I 


287 


11 


0.17 


35 03 


3 71 


-6.9( -10.7) 










4U1907+097 




HMXB.PSR 


C(B) 


287 


41 


9.83 


43. 74 


0.48 


10. 0{ 51.1) 


8.6 ± 1.1 


5.1 ± 0.9 


2.8 zb 1.1 


1.4 zb 2.8 


AQLX-l 


AQL X-1 


LMXBjNS 


At 


287 


82 


0.58 


35,72 


-4.14 


22. 8( 36.6) 


15.9 zb 1.1 


13.4 it 0.9 


9.3 zb 1.1 


1.8 zb 2.6 


SS433 




HMXBjNS? 




287 


96 


4.98 


39.70 




21. 7( 39.9) 


17.7 i 1.0 


11.8 i 1.0 


4.9 ± 1.0 


0.0 i 2.8 




Gratiat 1915+105 


TiTVTX Ft BHC 


At 








45.40 


-0.23 




1 -1- 1 1 
XOO.O ^ J. • X 






33.5 i 2.7 


4U1916-053 




LMXB,]MS 


B(N) 


289 


70 


-5.24 


31.36 


-8.46 


12. 5( 0.4) 


( 0.9 ± 0.7) 


( -0.2 ± 0.6) 


( -1.0 ± 0.6) 


( -1.6 ± 1.7) 


1940-04 


XB 1940-04 


LMXB.NS 


N 


295 


65 


-3.88 


35.32 


-13.15 


-1.8( -2.5) 


-1.5 ± 0.8 


0.6 ± 0.8 


-0.3 ± 0.8 


0.7 ± 2.2 


n n qa.aA-'xk 
Vj XV J J. 




HMXB PSR 


At 






^nnn 


63.20 


1.40 








" d's % 








n ivl 2s^ D J o rv 




9QS 


nn 


^9 ss 


66.19 


1.53 


"i* 2( ^iV) 


'^0^1 ± "I'l 


''o's ± *i" 1 


" ^ n r\ 








PSR 


X 


9ns 






68.77 


2,82 


9 Q ^ 9 \ 
- J. .Ciy z . o y 


9 Q -1- n t \ 


1 C\ f-\ -X- {'\ '7"\ 

U,u zt \i . i J 






H1954-|-319 




H M X B , N S 


/ 






; ■ 


68.32 


1,88 








1 f"*^ ±11 


-2 




CYG X 1 


srx. iVl j*v ±j , ji3 Jri 


" A ' 






■ 


71.34 


3.07 




fiT*^ Q -1- n Q 
oro.y zt u.y 


TniQ PI -1- HQ 

j.u±y.o zt u.y 


i-X^o.tzi zt u.y 


iJ^4.0 Zt i.O 


A. O iVl *i u U U -|- O 


onnn_i_oc; 

ij-inga 


J_j IVi D . JTL 










63-20 


-2.91 


0.1( -2.2) 


0.0 ± 1.1 


0.0 zb 0.9 


0.2 zb 1.1 


0.0 zb 2.8 


X T E J 2 1 2 -|- 38 1 




LMXB ,BHC 


1(C) 








75.39 


2.25 


-3.S( 3.2) 










W rtl4U 




n ivl A. j3 , IN o 




305 


12 


43 85 


80.93 


4.18 




1 9 ± 1 


3 6 ±0 9 


6 6 ±0 9 


11 3 ± 2 5 


GS2023+338 


Ginga 2023+338 


LMXB, BHC 


I(N) 


306 


02 


33.87 


73.12 


-2.09 


-7.5( 2.1) 


■ . . . ■ 


' . . . ■ 


' . . . " 


■ . . . 


NOVA.CYG 


Nova Cyg 1992 


Nova 


N 


307 


63 


52.63 


89.13 


7.82 


2.0( 2.9) 


1.0 ± 1.0 


1.1 ± 0.9 


-0.1 ± 0.9 


2.6 ± 2.2 


EXO2030+375 




HMXBjPSR 


At 


308 


06 


37.64 


77.15 


-1.24 


-7.0( -2.4) 


( 1.7 ± 0.7) 


( -1.0 ± 0.7) 


( -3.6 ± 0.7) 


( -5.3 zb 1.9) 


CYGX-3 


CYG X-3 


HMXBjBH? 


A 


308 


11 


40.96 


79.85 


0.70 


197. 8( 264.3) 


131.7 ± 0.8 


77.0 ± 0.8 


31.8 ± 0.8 


15.5 zb 2.1 


AEAQR 


AE AQR 


CV 


N(I) 


310 


04 


-0.87 


45.28 


-24.42 


-2.3( -3.2) 


-1.4 ± 0.8 


-0.6 ± 0.7 


-2.1 ± 0.7 


0.7 ± 2.0 


GRS2037-404 


Granat 2037-404 


LMXB, BHC? 


N 


310 


11 


-40.24 


1.46 


-37.28 


-0.6( -0.9) 


-0.6 ± 0.8 


-0.1 ± 0.7 


0.3 ± 0.7 


2.5 ± 2.1 


GROJ2058+42 




HMXB?,PSR 


At 


314 


23 


42.70 


84.04 


-1.79 


7.9( 2.0) 


( 1.7 ± 0.7) 


( 1.2 ± 0.6) 


( -1.4 ± 0.6) 


( -0.5 ± 1.7) 


SAXJ2103+45 


SAX J2103. 5+4545 


HMXB.PSR 


C 


315 


89 


45.75 


87.13 


-0.68 


6.8( 6.1) 


4.2 ± 0.9 


2.6 ± 0.7 


-0.4 ± 0.7 


3.4 ± 2.0 


4U2127-I-119 




LMXB,NS 


B 


322 


49 


12.17 


65.01 


-27.31 


11. 3( 15.8) 


5.6 ± 0.7 


5.2 ± 0.7 


2.8 ± 0.7 


-3.3 ± 2.0 


4U2129+470 


4U 2129+47 


LMXB,NS 


N 


322 


86 


47.29 


91.58 


-3.04 


0.8( 2.4) 


0.7 ± 0.9 


0.0 ± 0.7 


0.8 ± 0.7 


0.3 ± 2.0 


CEPX-4 


CEP X-4 


HMXB,PSR 


N(C) 


324 


90 


57.00 


99.03 


3.32 


1.1( 4-8) 


0,7 ± 0.9 


0.4 ± 0.7 


0.5 ± 0.7 


-0.7 zb 2.1 


SSCYG 


SS CYG 


DK 


C(B) 




68 


43.59 


90.56 


-7,11 


9.3( 13.2) 


5,7 ± 0,7 


3.4 zb 0.6 


3.0 zb 0.7 


1.0 zb 1.9 


CYGX-2 


CYG X-2 


L]VIXB,NS 


B 


326 


17 


38.32 


87.33 


-11.32 


19. 8( 27.7) 


15.0 ± 0.8 


3.7 zb 0.7 


2.6 zb 0.7 


2.0 zb 2.0 


PKS2155-30 


QSO 2155-304 


BL Lac 


C 


329 


72 


-30.23 


17.73 


-52.25 


6,3( 8,8) 


3,3 ± 0,8 


2.5 zb 0.7 


2.8 zb 0.7 


2.6 zb 2.0 


ARLACERT 


AR Lac 


RS CVn 


N 


332 


17 


45.74 


95.56 


-8.30 


-1,4( -1.1) 


-1,0 ± 0,7 


0.0 zb 0.7 


-0.2 zb 0.7 


2.6 zb 1.8 


G109-1 


GAL 109.1-01.0 


SNR,PSR 


N 


345 


28 


58.88 


109.09 


-1.00 


0.4( 2.2) 


0.1 ± 0.9 


0.3 ± 0.7 


0.5 ± 0.7 


0.9 ± 2.1 


CAS_A 




SNR 


C(B) 


350 


85 


58.81 


111.74 


-2.13 


7.2( 11.0) 


4.2 ± 0.9 


2.9 ± 0.7 


1.9 ± 0.7 


2.6 ± 2.0 



^Daily 4-band light curves in Crab units for the 83 category A & B sources and FITS files containing daily 20-100 keV fluxes and 16-channel spectra for all 179 sources are available on 

http: / / gam mar ay. msfc. nasa.gov/batsc /occult at ion/. 

denotes a category A transient source, otherwise the source is considered a persistent source. 

'^isted positions should be used for reference only. These positions were used in the analysis, and many are based On historical measurements and are often Only accurate to within a few arcminutes. Location 
errors of a few arcminutes do not significantly affect flux measurements with the BATSE EOT. 

^Conversion factors from mCrab units to photon and energy fluxes are given in Appendix B, Table B2. 

^Parenthesis denote cases where the uncorrected flux was very small relative to the systematic flux correction. In these cases the uncorrected 4-band fluxes are reported because we believe they are closer to the 
true 9.1-year averages. See section 4.3 for details. 



Table 5. BATSE Deep Sample Category "B" Sources 



Name 


Type 


Behavior 


F 


cor 








(ph cm 


-^s-i) 


SMCX-1 


HMXB,PSR 


P 


4.27E-03 


O.OOE+00 


NGC1275 


SY 




1.84E-03 


O.OOE+00 


XPER 


HMXB,PSR 




3.90E-03 


O.OOE+00 


LMCX-4 


HMXB,PSR 


P 


5.60E-03 


O.OOE+00 


4U0614+091 


LMXB,NS 


P 


5.98E-03 


-5.87E-04 


EXO0748-676 


LMXB,NS 




3.16E-03 


O.OOE+OO 


VELA 


SNR.PSR 




3.47E-03 


-2.57E-04 


1E1145-614 


HMXB,PSR 




6.06E-03 


8.64E-05 


3C273 


QSR 


P 


4.13E-03 


O.OOE+OO 


4U1323-619 


LMXB,NS 




4.00E-03 


1.08E-03 


NGC5548 


SYl 




2.22E-03 


O.OOE+OO 


PSR1509-58 


SNR,PSR 




3.76E-03 


1.21E-03 


4U1538-522 


HMXB,PSR 




3.66E-03 


-4.37E-04 


4U1627-673 


LMXB,PSR 




3.31E-03 


O.OOE+00 


4U1636-536 


LMXB,NS 




3.93E-03 


1.05E-03 


GX349+2 


LMXB,NS 




4.53E-03 


-6.12E-04 


4U1702-429 


LMXB,NS 


P 


4.31E-03 


-2.11E-03 


4U 1705-440 


LMXB,NS 


P 


5.51E-03 


-1.96E-03 


GX9+9 


LMXB,NS 




2.40E-03 


-1.47E-03 


GX354+0 


LMXB,NS 


P 


9.56E-03 


5.77E-04 


GX3+1 


LMXB,NS 




6.36E-03 


2.25E-03 


4U1812-12 


LMXB,NS 




8.13E-03 


-3.26E-04 


4U1820-303 


LMXB,NS 




1.55E-03 


-1.68E-03 


H1822-371 


LMXB,NS 




4.06E-03 


O.OOE+OO 


LS5039 


HMXB,BHC 




1.20E-03 


-1.43E-03 


SS433 


HMXB,NS? 




6.36E-03 


1.20E-03 


4U1916-053 


LMXB,NS 




6.27E-05 


-2.61E-03 


4U2127+119 


LMXB,NS 




1.87E-03 


O.OOE+OO 


CYGX-2 


LMXB,NS 




3.40E-03 


O.OOE+OO 



SF tear (20-40) (40-70) (70-160) (160-430) 



(ph 


cm 


VI) 








(mCrab) 




1 


06E- 


04 


1 


40 


24.2± 0.6 


9.2+ 0.5 


-0.6+ 0.6 


-.70+ 1.5 


1 


23E- 


04 


1 


40 


7.1+ 0.6 


3.1+ 0.5 


3.3+ 0.5 


3.9+ 1.4 


1 


17E- 


04 


1 


40 


9.6± 0.5 


12.7± 0.5 


6.5± 0.5 


-0.2+ 1.3 


1 


55E- 


04 


1 


40 


20.2± 0.6 


9.4± 0.6 


0.7± 0.6 


-8.0+ 1.5 


1 


llE- 


04 


1 


47 


20.6+ 0.6 


16.6+ 0.5 


12.4+ 0.5 


8.0+ 1.4 


1 


46E- 


04 


1 


40 


9.2+ 0.7 


8.8+ 0.6 


5.1+ 0.6 


-6.1+ 1.6 


1 


42E- 


04 


1 


77 


9.5+ 0.6 


9.3+ 0.6 


10.6+ 0.6 


11.6+ 1.6 


1 


62E- 


04 


1 


62 


20.5+ 0.7 


14.8+ 0.6 


4.2+ 0.7 


-6.5+ 1.8 


5 


71E- 


05 


1 


40 


18.4+ 0.5 


19.2+ 0.5 


23.1+ 0.5 


32.3+ 1.4 


1 


46E- 


04 


1 


66 


11.3+ 0.7 


11.2+ 0.6 


6.3+ 0.6 


-0.3+ 1.6 


1 


19E- 


04 


1 


40 


6.3± 0.5 


6.3± 0.5 


5.7± 0.5 


5.8+ 1.4 


1 


42E- 


04 


1 


53 


10.3+ 0.6 


10.3+ 0.6 


11.7+ 0.6 


14.1+ 1.6 


1 


49E- 


04 


1 


50 


14.8+ 0.7 


6.8+ 0.6 


1.1+ 0.6 


-3.8+ 1.7 


1 


60E- 


04 


1 


40 


14.5+ 0.7 


4.2± 0.6 


-1.3± 0.7 


-7.3± 1.8 


1 


84E- 


04 


1 


44 


15.8+ 0.8 


7.4+ 0.7 


0.7+ 0.8 


-5.1+ 2.1 


2 


65E- 


04 


1 


35 


20.3+ 1.2 


5.0+ 1.1 


2.6+ 1.1 


0.0+ 3.0 


4 


02E- 


04 


1 


30 


14.1+ 1.8 


11. 2± 1.6 


10.0+ 1.7 


9.4± 4.7 


2 


78E- 


04 


1 


31 


17.6+ 1.2 


10.3+ 1.0 


7.9+ 1.1 


12.7+ 2.9 


2 


14E- 


04 


1 


49 


9.2+ 1.0 


3.8+ 0.8 


2.2+ 0.9 


5.2+ 2.3 


1 


40E- 


04 


1 


28 


49.4± 1.0 


37.1± 0.9 


19.1+ 0.9 


1.9+ 2.4 


2 


20E- 


04 


1 


43 


23.4+ 1.0 


12.4+ 0.9 


8.3+ 0.9 


12.9+ 2.4 


3 


56E- 


04 


1 


47 


20.7+ 1.6 


23.5+ 1.4 


20.4+ 1.5 


14.3+ 4.0 


1 


70E- 


04 


1 


37 


10.6± 0.8 


-1.7± 0.7 


-3.4± 0.7 


-8.5+ 1.8 


1 


64E- 


04 


1 


40 


21.2+ 0.7 


2.5+ 0.6 


-2.2+ 0.7 


-6.5+ 1.8 


1 


66E- 


04 


1 


48 


3.2+ 0.7 


2.9+ 0.7 


2.0+ 0.7 


3.7+ 1.8 


1 


59E- 


04 


1 


49 


21.9± 0.7 


14.7± 0.6 


6.1+ 0.7 


-0.2+ 1.8 


1 


48E- 


04 


1 


44 


0.6+ 0.7 


-0.4+ 0.6 


-1.0+ 0.6 


-1.6+ 1.7 


1 


18E- 


04 


1 


40 


5.7+ 0.5 


5.1+ 0.5 


2.8+ 0.5 


-2.6+ 1.4 


1 


23E- 


04 


1 


40 


15.2± 0.6 


3.8± 0.5 


2.5+ 0.5 


2.0+ 1.4 



Note. — Definitions for F, 5F, Fcor, and ecor are given in Eq. (1). 



I 

CO 

-J 
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Tablc 6. Orbital. Siipcrorbital and Quasi- Periods Observed in BATSE EO Light Cnrves 



Other Instruments/Meth Earth Occultation 



Name 


Porb (d) 


Plong (d) 


refs 


FFT Peaks 


FAP 


Category A 


4U 1700-377 


3.411581 




[4] 


3.411±0.001 


3.43E-65 


CEN X-3 


2.08706533 


~140 


[1] 


2.0873±0.0003 


5.28E-25 


CYG X-1 


5.599829 


294 


[3,8] 


5.6" 




CYG X-3 


4.792426h 




[6] 


4.79249h'' 




EXO2030+375 


46.0 




[7] 


45.9±0.1 


5.78E-20 


GRO J2058-H42 


~110 




[2] 


97.5±0.4 


2.88E-3 


GRS 0834-430 


105.8 




[2] 


broad power 97-118 


<7.48E-6 


GS 1845-024 


242.18 




[5] 


241.9±4.5 


1.63E-10 


GX 301-2 


41.498 




[2,3] 


41.49±0.07 


4.91E-77 


HER X-1 


1.700167412 


35 


[1,2] 


34.80±0.04 


<1.0E-99 


OAO 1657-415 


10.44809 




[2] 


10.45±0.01 


2.21E-46 


VELA X-1 


8.964368 




[2] 


8.963±0.006 


<1.0E-99 


Category B 


4U 1538-522 


3.72839 




[2] 


3.728±0.001 


1.43E-2 


LMC X-4 


1.40841 


30.4 


[1] 


30.34±0.03 


5.63E-52 


SMC X-1 


3.892116 


-55 


[1,3] 


3.892±0.001 58.6±0.2 


9.95E-2,3.37E-7 



References. — [1] Wijers & Pringle (1999) [2] Bildsten et al. (1997) [3] Priedhorsky & Holt (1987) [4] Rubin 
et al. (1996) [5] Finger et al. (1999) [6] Singh et al. (2002) [7] Wilson et al. (2002) [8] Brocksopp et al. (1999a) 



''Detected only in FFTs of single step data 
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Table 7. Earth Occultation Imaging Results for Category B Sources 



Source Name 


# of Days 


Channel/Position Results 






visible ill all 3 channels, some positions 




1 A 


visible in all 3 channels, some positions 




16 


slight indication of excess in some positions 


IjMU A- 4 


10 


strong in ch 1-2, visible in 3, most positions 


4U UDi4+Uyi 


1 Q 

lo 


visible in ch 1 when in center of grid 




1 Q 

13 


strong in ch 1-2, most positions, sometimes in ch 3 


Vela pulsar 


11 


no convincing image — noisy sky around source 


IT? II/IK C-l A 

lH/ ii4o-Oi4 


1 /I 

14 


slight indication of excess in some positions, ch 1 




11 


visible in some positions, strongest in ch 3 


/ITT 1 QOQ (^1 O 

4U iozo-Diy 


lo 


strong in most positions, chs 1-3 


NGC 5548 


13 


no convincing image 


PSR 1508-59 


15 


visible in all 3 channels, some positions 


4U 1538-422 


12 


visible only in center of grid 


4U 1626-673 


10 


visible in some positions, ch 1 only 


4U 1636-536 


10 


visible in some positions, ch 1-3 


GX 349+2 


13 


visible in most positions, ch 1-2 


4U 1702-429 


11 


confusion with 4U 1705-440, ch 1-3 


4U 1705-440 


12 


strongly visible most positions - confusion with 4U 1702-429 


GX 9+9 


15 


no convincing image 


GX 354+0 


14 


visible most positions, ch 1-3 


GX 3+1 


14 


visible in center ch 1-3, else sky too noisy 


4U 1812-12 


11 


bright, distinguishable from GX 17+2 in ch 2 


4U 1820-303 


14 


no convincing image 


H 1822-371 


15 


visible some positions, ch 1-2, not compact 


LS 5039 


12 


slight excess in some positions, ch 1-3 


SS 433 


12 


slight excess in some positions, ch 1 


4U 1916-053 


13 


no convincing image 


4U 2127+119 


14 


no convincing image 


CYG X-2 


14 


no convincing image 
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A. Special Notes on Transient Activity in 
the BATSE Database 

The available literature on the results of BATSE 
detection and monitoring for transient sources is not 
always complete. The difficulty of identifying new 
sources, in particular, directly stems from the lim- 
ited spatial resolution and sensitivity of the BATSE 
EOT. As discussed in Harmon et al. (2002), our abil- 
ity to spatially resolve point sources on the sky is 
limited by the geometry of the Earth's limb at the 
time of measurement, coupled with the slow change 
of this geometry (~days) as the spacecraft orbit pre- 
cesses. This is compounded by the sometimes brief 
duration of transient outbursts. Some special notes 
are provided here regarding activity that was not as 
well understood at the time of its occurrence. 

EXO 1846-031, GS 1843+00 and GS 1845-024 

These sources, occupying a relatively small (few- 
degree) region on the sky, were all detectable at var- 
ious times in the BATSE database. First noted in 
a deep search of the BATSE database using epoch 
folding, a source with a 241 day period was given the 
name GRO J1849-03 (Zhang et al. 1996a, 1996b). At 
that time, the regularity of the outbursts suggested 
that GRO J1849-03 was a HMXB, but pulsations 
were not detected. Later sensitivity enhancements in 
BATSE pulsar analysis, however, revealed 94-s pul- 
sations detected at times of outbursts. This obser- 
vation, combined with a pulse arrival timing analysis 
that revealed an orbital period of 242. 18 ± 0.01 days, 
identified GRO J1849-03 with GS 1845-024 (Finger 
et al. 1999). A more precise positional confirmation 
was also made with the Wide Field Camera on Bep- 
poSax (Soffitta et al. 1998). 

BATSE observed short-lived transients with a 
hard spectrum in 1993 (Harmon et al. 1993) and 
slightly less than a year later (Zhang et al. 1994) 
in this same region of the sky. The spectral shape, 
both times, was reminiscent of a BHC. Unfortunately, 
the limb geometry was not sufficiently to locate the 
source to more than a few degrees in one direction 
within the few days that the outbursts were observ- 
able. A follow-up observation two weeks after the 



1994 outburst did not indicate any unusual optical 
activity in the vicinity of EXO 1846-031 (Grindlay, 
Garcia, & Zhao 1994). The BATSE error boxes were 
consistent with the reported discovery location of the 
black hole transient EXO 1846-031 (Parmar et al. 1993), 
and so we tentatively identified the outbursts with 
this source (Zhang et al. 1994). 

In Table 3, we include the regular outbursts of 
GS 1845-024, a single outburst of the pulsar GS 1843+00 

(Wilson et al. 1997b) and the two outbursts with hard 
spectra attributed to EXO 1846-031. 

Galactic Center Transient GRS 1739-278 

A transient source was first detected on 16 Mar 
1996 with Granat/Sigma. (Paul et al. 1996) in the 
Galactic center region with a hard spectrum. A study 
of its properties with Granat (Vargas et al. 1997) in- 
dicated that the source, named GRS 1739-278, prob- 
ably contained a black hole. Because of its proximity 
(-^1°) to a bright source discovered with BATSE only 
a few months earlier (GRO J1744-28 or "bursting pul- 
sar") (Fishman et al. 1995; Kouveliotou et al. 1996), 
it was difficult to separate the fluxes from the two 
sources. The first observations showing a clear detec- 
tion of GRS 1739-278 by RXTE/ASM was on TJD 
10143 (1 Mar 1996), although a single measurement 
exists on TJD 10101 (19 January 1996) that indicates 
about a 3a detection. Greiner, Dennerl, & Predehl 
(1996) made a follow-up observation to the Granat 
discovery on TJD 10158 (18 Mar 1996) to get a more 
accurate location for the X-ray counterpart, and dis- 
covered a substantial X-ray scattering halo. Due to 
the delay in arrival time of the scattered photons, the 
ROSAT results implied that the source could have 
first appeared as early as 1995 November-December. 
A very sparsely sampled BATSE light curve for GRS 
1739-278 indicated the source was visible as early as 
~10035 (14 Nov 1995) at a level of ~200 mCrab in 
the 20-100 keV band, and decaying below the source 
confusion limit (^50 mCrab) in the Galactic center 
region by ~10497 (18 Feb 1997). These beginning 
and ending dates for the outburst were entered into 
the flare database and Table 3. 

RXTE Transients XTE J1550-564, XTE J1859-I-226, 
XTE J1118-h480, XTE J1543-568, XTE 30111-7317 
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and V464I SGR (Sagittarius) 

These transients were discovered in outburst by 
RXTE during the final years of the CGRO mission. 
Several of these transients were quite bright and 
reached a Crab or more in the BATSE energy band- 
pass. The RXTE/ ASM, which was designed to quickly 
locate and monitor new transients, generally was able 
to report these outbursts rapidly, sometimes even 
before BATSE data containing the same transient 
were received for processing. It is interesting scientifi- 
cally, as well as important to our accounting of bright 
sources, to examine archival data for these transients, 
especially prior to the launch of RXTE in 1995. 

As part of the deep sample, light curves were 
generated for XTE J1550-564, XTE J1859+226, XTE 
J1118-I-480, XTE J1543-568, XTE JOlll-7317 and 
V4641 SGR. All showed detectable outbursts as re- 
corded in Table 3. Examination of light curves from 
the beginning of the CGRO mission in April 1991 
prior to the known outbursts occurring at various 
times after December 1995 did not show any emis- 
sion more than five successive days greater than ~100 
mCrab intensity in the 20-100 keV band. Outbursts 
of lesser duration and intensity could possibly be pre- 
sent in the dataset, but require more in-depth analy- 
sis. 



B. Energy Channel Mapping Scheme, Flux 
Conversion Factors cind Data Problems 

Pulse height data from the BATSE large area de- 
tectors were normally acquired into 128 high energy 
resolution (HER) channels. Duo to telemetry limits, 
the HER data were mapped into 16 broader energy 
channels and read out at 2.048s to create the contin- 
uous (CONT) background data. EOT measurements 
were made with the CONT dataset. Normally the 
CONT channel bin edges were kept to a constant set 
of values, that, coupled with the LAD gain control, 
kept the effective energy loss pulse height response 
per channel extremely stable. However, on rare occa- 
sions, the 128 to 16 channel mapping was changed to 
accommodate a greater energy dispersion when ob- 
servations of sources with softer spectra and possible 



line features were desired (e.g., accreting pulsars or 
soft gamma-ray repeaters). Table Bl shows several 
mapping schemes (numbered 0, 1, 2, 3, 4a and 4b) 
that were used during the mission and the times in 
which these were applied. From column 2, it can be 
seen that scheme 2 was used about 90% of the time. 

Schemes 2, 3, 4a and 4b had in common several 
of the same HER channels for which we could obtain 
4-energy band data spanning changes in the HER to 
CONT channel mapping. Using the Crab fluxes as a 
standard candle, a tool was created for which the flux 
of any source relative to the Crab could be generated 
in the energy bands 20-40, 40-70, 70-160 and 160-430 
keV. The Crab relative rates quoted in Tables 4 and 
5 can be converted to photon number per keV (ph 
cm~^s~^ keV~^) or energy flux per band using the 
multipliers given in Table B2. Note that this method 
of four-band conversion to Crab relative rates cannot 
be used for schemes or 1, and so data before TJD 
8406/81139 and for the brief intervals TJD 8807/8279 
- 8809/64510 and 8812/65817-8812/66102 cannot be 
treated in this manner. 

Table B3 is a list of days when unusual or anoma- 
lous events occurred during the CGRO mission and 
no Earth occultation measurements could be extracted. 
These events produced gaps in most or all source flux 
histories at the given times. Routine events causing 
loss of measurements such as source interference, ex- 
ceeding the critical angle for Earth occultation, or 
telemetry gaps are not included. Table B3 gives the 
times, data condition reported, and the cause of the 
anomaly. 

C. Comparison to JPL Enhanced BATSE 
Occultation Package (EBOP) 

Skelton et al. (1994), Ling et al. (1996, 2000) 
and co-workers at the Jet Propulsion Laboratory de- 
veloped an Earth occultation software package over 

a period of several years with the intent to generate 
higher sensitivity results than can be obtained with 
the treatment described in Harmon et al. (2002). The 
JPL package is known as the Enhanced BATSE Oc- 
cultation Package (EBOP) and has been applied to 
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BATSE data independent of the efforts at MSFC. 
Ling et al. (2000) published a compendium of mea- 
surements (also called a catalog) for the 1991 May 
to 1994 October epoch. We feel it is very impor- 
tant to compare the results from EBOP with those 
of this work, since there are significant discrepancies 
between what we report and what can be found in 
Ling et al. (2000). 

Primarily, there are two major differences be- 
tween the JPL method and the method used here. 
First, EBOP uses a semi-physical model for the detec- 
tor background counting rates. The model is based 
on expected contributions of low-energy gamma-ray 
fluxes local to the low Earth orbit (LEO) environ- 
ment. These include cosmic ray secondary radiation 
and activation products from orbital passes through 
radiation fields in LEO. The JPL global background 
model consists of a mix of the local radiation com- 
ponents and a combination (determined by the fit) 
of the BATSE spectroscopic detector (Fishman et 
al. 1989a, b) counting rates as a predictor of the low- 
energy background in the LADs. Secondly, the ex- 
traction of source signals is usually performed in one- 
day segments, with a single fit including terms for all 
sources in the EBOP catalog with no a priori assump- 
tions of their intensity. The main goal of this work is 
increase sensitivity by fitting much longer intervals of 
background data than the four minutes typically used 
in our work, which does not require a sophisticated 
background model. 

In Table CI, we compare selected results for a 
few sources from the deep sample in Table 4 and the 
EBOP catalog (Ling et al. 2000). To best illustrate 
the comparison, wc selected sources over a range of 
flux intensities, as well as to pick sources that did 
not show a large degree of variability on timescales of 
months to years. Table CI shows fluxes in three se- 
lected energy bands with centroid energies that over- 
lap in EBOP and the Table 4 deep sample. Results 
from the deep sample are quoted with and without 
systematic error corrections. The EBOP results arc 
reported as minimum and maximum average fluxes 
for for several ^400 day intervals in Table 3 of Ling et 
al. (2000). The deep sample results are the 9.1y aver- 
ages from this work (Table 4) converted from mCrab 



to photon flux with the multiplicative factors from 
Table B2. For the Crab, the two methods give very 
similar results, and confirm the cross-check of the of 
the two methods from Ling et al. (2000). Sources 
of moderate intensity, such as 4U 1700-37 and NGC 
4151, show significant discrepancies, particularly for 
the highest of the three energy bands. For relatively 
weak sources at the few mCrab level, such as the 
LMXB Circinus X-1, the supernova remnant PSR 
1509-58, and the Seyfert Galaxy NGC 1275, our av- 
erage fluxes are roughly an order of magnitude less 
than the EBOP results. For Set X-1, wc only show 
upper limits (<2(t) (corrected for systematic error), 
whereas the JPL method yielded significant broad- 
band emission at least an order of magnitude higher. 

To make sure that the time interval for which 
the averages were constructed was not a major factor 
in the comparison, we generated long term spectra 
for TJD 8393-8800 as was done in Ling et al. (2000). 
Figures 17a and 17b show the broad-band spectra 
of Cir X-1 and and PSR 1509-58 using the MSFC 
method with no systematic error correction applied, 
the EBOP results from Ling et al. (2000), and the Ta- 
ble 4-band average fluxes with systematic error cor- 
rections. The broad-band spectra from the JPL and 
MSFC methods are considerably different over most 
of the sensitive energy range of the LADs. Differences 
between the broad-band MSFC data and the 4-band 
averages for Cir X-1 are mostly caused by the large 
systematic error correction for that source. The find- 
ings reflect the same conclusions made using Table 
CI. 

Ling et al. (2000) discuss sources of error for the 
Earth occultation analysis specific to EBOP in some 
detail and point out that the high energy fiuxes they 
see for some sources such as Set X-1 and NGC 1275 
should be viewed with caution. Since we observe con- 
siderably lower fluxes for these and other sources in 
the 160-430 keV band, it suggests that their concern 
could be justified. 

Here we point out possible reasons for the dis- 
crepancy between our work and the EBOP results. 

One possibility concerns unaccounted for components 
in the gamma-ray background for which the LADs 
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havc significant sensitivity, but are not included in 
the JPL model. We have found that, when fitting 
the background near occultation steps for more than 
a few minutes, there is an increased likelihood for 
slow, very significant variations to cause the mea- 
sured background to be considerably different from 
the assumed model. These variations are not associ- 
ated in any way with occultation steps. Such effects 
could be high energy diffuse emission from the Galac- 
tic plane or unusual local background variations on 
specific days (e.g., solar activity or electron precipita- 
tion events) that can affect the computed level of the 
background in the high energy (a few hundred keV) 
data. These variations from the model can be trans- 
ferred to source terms which are fit simultaneously 
with the backgroimd. Another possibility is the as- 
sumed set of sources used by the JPL method. JPL 
used a set of 64 sources, all of which are fit simul- 
taneously, regardless of any intensity or time infor- 
mation. Our catalog consisted of 179 sources, how- 
ever, no more than 10 sources were fit simultaneously 
in the MSFC method, and usually only 1-3 per 4- 
minute fitting window were typical. This reduces the 
tendency toward large scale coupling between source 
terms, and presupposes that we have partial informa- 
tion about what sources a detector is likely to be see- 
ing. This may not explain anomalously large fluxes 
at high energy in the JPL data, but it can help to 
explain the tendency to yield fluxes higher than the 
MSFC method for lower intensity sources. Although 
a finite source set is also a problem for the MSFC 
method, the strategy was to iteratively improve char- 
acterization of the sources in the catalog based on 
average brightness, times of transient activity, and to 
use a source set taken from published catalogs. Our 
method therefore yields improved results by having a 
better knowledge of source activity each time a pass 
through the dataset is performed. 
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Table Al. LIST OF ACRONYMS AND ABBREVIATIONS 



Acronym or Abbrev. Meaning 



AGN active galaxy or galactic nucleus 

ASM All-Sky Monitor 

BATSE Burst and Transient Source Experiment 

BH or BHC black hole candidate 

BL Lac BL Lacertae object 

CGRO Compton Gamma Ray Observatory 

CH channel 

CONT large area detector continuous data 

CV cataclysmic variable 

DN dwarf nova 

EBOP JPL Enhanced BATSE Occultation Package 

EGRET Energetic gamma-ray Experiment Telescope 

EXIST Energetic X-Ray Imaging Survey Telescope 

EOT Earth occultation technique 

FFT fast Fourier transform 

GCR Galactic center region 

HEAO High Energy Astronomy Observatory 

HER high energy resolution 

HISGEN BATSE photon flux history generator 

HMXB high mass x-ray binary 

INTEGRAL International Gamma-Ray Astrophysics Laboratory 

JPL Jet Propulsion Laboratory 

LAD BATSE large area detector 

LMXB low mass x-ray binary 

MSEC Marshall Space Flight Center 

NS neutron star 

OCCSYS automated data flagging software 

OF_W MSFC EOT flux history generator 

OTTB optically thin thermal bremsstrahlung spectral fitting model 

PL power law spectral fitting model 

PSR pulsar 

QSR quasar 

ROSAT German X-Ray Roentgen Satellite 

RS CVn eruptive variable of RS Canum Venaticorum type. 

RXTE Rossi X-Ray Timing Explorer 

S/C spacecraft 

SGR soft gamma-ray repeater 

SNR supernova remnant 



Tabk^ Al ContiiiiKxl 



Acronym or Abbrev. Meaning 



SOD seconds of day 

SOI source of interest 

SY(1 or 2) Soyfcrt galaxy, type 1 or 2 

TJD Truncated Julian Date (Julian Date - 2,440,000.5) 

UVS ultraviolet source 

XRS X-ray source 



Table Bl. BATSE LAD 128-Channel to 16-Channel Mapping (Lookup Table) 



Scheme # 


%Time 


TJD 

(Days) 


SOD 

(sees) 





1 


2 


3 


4 


5 


CONT energy 
6 7 8 


channel 
9 


10 


11 


12 


13 


14 


15 





0.2 


8361 








13 


15 


17 


19 


22 


25 


30 


37 


49 


64 


74 


89 


100 


111 


125 


1 


1.2 


8367 








9 


11 


14 


18 


22 


26 


30 


38 


49 


65 


75 


90 


98 


106 


122 


2 


12.0 


8406 


81139 





7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 





0.1 


8807 


8279 





13 


15 


17 


19 


22 


25 


30 


37 


49 


64 


74 


89 


100 


111 


125 


2 


0.1 


8809 


645 1 1 


Q 


7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 





0.0 


8812 


65817 





13 


15 


17 


19 


22 


25 


30 


37 


49 


64 


74 


89 


100 


111 


125 


2 


17.7 


8812 


66103 





7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 


3 


0.1 


9400 


59319 


7 


9 


11 


14 


18 


20 


23 


25 


28 


32 


36 


49 


65 


75 


90 


122 


2 


19.3 


9419 


60041 





7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 


4a 


3.5 


10062 


71275 


122 


7 


8 


9 


10 


11 


12 


14 


16 


18 


23 


28 


36 


49 


75 


90 


2 


8.6 


10178 


85302 





7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 


4b 


1.2 


10465 


83335 





7 


8 


9 


10 


11 


12 


14 


16 


18 


23 


28 


36 


49 


75 


122 


2 


14.4 


10504 


80991 





7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 


4b 


2.1 


10983 


71907 





7 


8 


9 


10 


11 


12 


14 


16 


18 


23 


28 


36 


49 


75 


122 


2 


0.1 


11053 


51166 





7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 


4b 


0.1 


11055 


54394 





7 


8 


9 


10 


11 


12 


14 


16 


18 


23 


28 


36 


49 


75 


122 


2 


18.2 


11085 


74663 





7 


9 


11 


14 


18 


23 


28 


36 


49 


65 


75 


90 


98 


106 


122 
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iable Bz 


Conversion Factors for 4-Energy Band Data in 


iabie 4 


HER Channels 


Energy Band 
(keV) 


Log Avg 
(keV) 


Multiplier 
(ph cm ~^s~^keV~^) 


Multiplier 
(ergs cm~^s~^) 


7- 11 

11 - 18 
18 - 36 
36- 75 


20-40 

40-70 

70-160 

160-430 


28.3 
52.9 
105.8 
262.3 


8.508 X 10"*^ 
2.323 X 10"^ 
5.514 X 10""^ 
6.261 X 10"^ 


7.567 X 10-12 
5.830 X 10-12 
8.079 X 10-12 
6.887 X 10-12 



Note. — ower HER channel is the lower energy edge of the corresponding CONT 
channel and the higher HER channel is the corresponding upper edge of the correspond- 
ing CONT channel as given in Table Bl. Example: HER channels 7 - 11, corresponds 
to CONT channels 1 - 2 in scheme 2 and CONT channels 1 - 4 in scheme 4a. 



Tabic B3. List of Mission Anomalies/Data Loss Events Affecting BATSE Earth Occultation Results 



TJD/SoD 


Condition 


Cause of Anomaly 


8449/05500-12000 


some LAD gains out of balance 


SAA passage w/HV on 


8490/03000-35000 


some LAD gains out of balance 


SAA passage w/HV on 


8588 


observing plan read error 


? 


8741/80355 - 8742/11519 


data crash 


flight software problem 


8742/63376 - 8742/83566 


data crash 


flight software problem 


8807/03065-66377, 68578-71280 


HV off 


S/C power supply problem 


8807/66300 - 8809/53160 


LAD gains out of balance 


S/C power supply (cont.) 


8812/65733-65964, 68510-71490 


HV off 


change S/C power supply 


8812/65960-68512 


LAD gains out of balance 


S/C power supply (cont.) 


9075/57503-62797 


data crash 


flight software problem 


9075/79431 - 9076/09796 


data crash 


flight software problem 


9110/84820 - 9112/54665 


HV off 


reboost activities 


9113/81403 - 9115/53135 


HV off 


reboost activities 


9152/84296 - 9154/16964 


HV off 


reboost activities 


9336/49077-65666 


HV off 


reboost activities 


9336/65666-75000 


some LAD gains out of balance 


reboost activities 


10014/03349-15300 


LAD gains out of balance 


SAA passage w/HV on 


10014/15300-65700 


LAD 3 gain out of balance 


SAA passage w/HV on; LAD 3 reset 


10014/15300 - 10016/10008 


LAD 6 gain out of balance 


SAA passage w/HV on; HV reset 


10416/16897-19200 


bad rates 


recovery from flight software upset 


10417/70332-71000 


bad rates 


recovery from flight software upset 


10416/16897 - 10418/69320 


LAD gains slightly out of balance 


flight software upset 


10900 - 11690 


LAD 7 gain decreased 


AGO caused HV to reach maximum value 


11001/11000-50000 


LAD gains out of balance 


SAA passage w/HV on 


11001/50000 - 11002/30400 


LADs 0,6 gains out of balance 


SAA passage w/HV on 


11002/30400 - 11003/46570 


LAD gain out of balance 


SAA passage w/HV on 


11003/46570-50500 


LADs 0,2,6 gains out of balance 


HV reset 


11046/74756 - 11047/28002 


data lost 


bad solar ephemcris 


11047/30897-71816 


data lost 


bad solar ephemcris 


11047/71816 - 11050/77609 


HV off 


bad solar ephemcris 


11307/04563-31857 


data lost 


high-gain antenna problems 


11354/03716 - 11359/72432 


data lost 


high-gain antenna problems 


11518/55813 - 11522/00089 


data lost 


gyro failure 



Table CI. Comparison of JPL EBOP and MSFC Catalog Results 



Source 


Energy Band 
(keV) 


EBOP 
ph cm~^s~^keV~^ 


Energy Band 
(keV) 


MSFC 
ph cm~^s~^keV~^ 


MSFC/w cor. 
ph cm~^s~^keV~^ 


CRAB 


45-55 


Z, .OO — O.J-UA ±U 


40-70 


o "^R V 1 

.i.OU A ±yj 


9 V 1 

. oo A lyj 




98-123 


'^ 10— '^ ROvIO"^ 

O . i-\J — O .\J\J i-\J 


70-160 


c, fii V 1 n~'^ 

O - vJ -L A ± U 


5 64 X 1()~^ 




.iOU-O xo 


f\ ^ A 7 ni V 1 

U.-L4: — I .U-LA ±U 


1 fin A'^n 


6 47 X 10~^ 


fi f^n V 1 n~ "'^ 
D . ou X ±u 


ft U -L ( UU-O ( 


45-55 


41. oy — o.oyx lu 


40-70 


"i AQ V 1 n"** 


"X A9i V 1 n~'^ 




98-123 


QQ A 'iA V 1 


70-160 


'i 1^9 y 1 


"x f^i vin~"'^ 

O . O -L A ±U 




230-313 


1 7n A AA V 1 


1 fin 4*^0 
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FIGURE CAPTIONS 

Fig. 1.— Comparison of sensitivity (5ct) in photons cm-^s-^keV-^ for BATSE, the 1977-1979 HEAO 1 A-4 
sky survey (Levine et al. 1984) and the survey concept mission EXIST (Grindlay et al. 2001, 2003). Nominal 
energy bins for the two historical instruments are shown. An approximate representation of the Crab Nebula 
spectrum is shown for reference. 

Fig. 2. — Aitoff-Hammer projection of the sky in Galactic coordinates showing the locations of the 179 
sources in the BATSE Earth occultation catalog. The sources arc labeled by object type. 

Fig. 3. — (a) The Compton Gamma Ray Observatory showing placement of the BATSE detector modules. 
The detector modules on the top four corners of the spacecraft are numbered 0,2,4, and 6, proceeding 
clockwise from the detector facing the reader above the high gain antenna. Similarly the bottom four 
modules are numbered 1,3,5, and 7.(b) Major components of the BATSE detector modules (one of eight). 

Fig. 4. — Flow diagram for processing of BATSE Earth Occultation data. Part (a) shows major processing 

steps beginning with components of the BATSE data archive on the far left (see text) leading to raw (counts 
per detector) rate histories. The raw histories were then converted to photon fluxes (20-100 keV) using the 
BATSE detector response, and also to intensities relative to the Crab Nebula flux in four flux bands (20-40, 
40-70, 70-160 and 160-430 keV). Part (b) shows the final processing steps (epoch folding, FFT analysis, 
deriving average flux and uncertainty) in generating the tables in the catalog. The "Crab model" allows 
conversion of raw count rates directly into Crab relative fluxes (see text). Iteration #n refers to a single pass 
through the entire 9.1y dataset, which can be repeated to increase the size of the sky sample and to reduce 
systematic error. 

Fig. 5.— Eclipse profiles in the 20-40 keV band for the five bright persistent HMXBs 4U 1700-377, Gen X-3, 

Her X-1, OAO 1657-415 and Vela X-1 in the BATSE dataset. Published ephcmcrides are used to eliminate 
fitting terms for these binaries while eclipses are in progress. Table 6 indicates the orbital periods used in 
epoch folding of the occultation flux histories shown in the flgure. 

Fig. 6. — Flux histories (20-100 keV, 1991 April - 2000 June) for four persistent black hole binary systems 
observed with BATSE. 

Fig. 7. — Flux histories (20-100 keV, various times) for flve transient black hole binary systems. 

Fig. 8. — Flux histories (20-100 keV, 1991 April - 2000 June) for five low mass X-ray binary neutron star 
systems. 

Fig. 9. — (a) Flux histories (20-100 keV, 1991 April - 2000 June) for six persistent accreting X-ray pulsars. 
For purposes of this catalog, a persistent source is defined as consistently exceeding 0.01 photons cm~^ s~^ 
(20-100 keV, ~ 35 mCrab). Two of these systems GX 1-1-4 (P^pin =i 120 s) and Her X-1 (P^pin 1.24 s) have 

a low-mass companion, while the rest, Cen X-3 (P^pin ^ 4.8 s), GX 301-2 (Pspin ^ 681 s), OAO 1657-415 
(Pgpin — 10.4 s), and Vela X-1 (Pgpin — 283 s) have a high-mass evolved supergiant companion, (b) Flux 
histories (20-100 keV, various times) for four transient X-ray pulsars. For purposes of this catalog, a transient 
source is defined as exceeding 0.01 photons cm^^ s~^ (20-100 keV, ~ 35 mCrab) for an identifiable, but 
limited period of time. All four sources shown, EXO 2030-^375 (Pspin ^ 42 s), GRO J1944+26 (Pspin ^ 15.8 
s), GRS 0834-430 (Pspin ^ 12.3 s), and GS 1845-02 (Pspi„ ~ 94.8 s), have an identified Be star companion. 



Fig. 10. — Flux histories (20-100 keV, 1991 April - 2000 June) for four galaxies with gamma ray active 
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nuclei. 

Fig. 11. — Sky grid, consisting of 162 points, used for determination of location-dependent systematic error. 
Grid spacing is described in the text. 

Fig. 12. — (Top) Measured averages in photons cm~^s~^ (20-100 keV) and (Bottom) width (in as or standard 
deviations) of the 9.1y flux histories for the grid locations shown in Fig. 11 as a function of Galactic longitude. 

Fig. 13. — Distribution of 9.1y flux averages (20-100 kcV) in order of brightness before (filled circles) and 
after correction with systematic error model. Symbols denote source categories used for the deep sample 
results. 

Fig. 14. - Histogram of fractional changes in flux significance of the deep sample when systematic errors 
were applied. The peak in the histogram at ~-0.3 corresponds to cases where the correction to the flux is 
zero, and the correction is dominated by the correction to the flux uncertainty. Systematic corrections to 
the average fluxes in Table 4 were only applied when the computed fractional change was between the two 
vertical dashed lines. 

Fig. 15. — (a) Flux histories in the indicated energy band for the three category B sources (LMC X-1, SMC 
X-1 and 4U 1538-522) detected in the FFT analysis and epoch folded at the appropriate orbital periods as 
given in Table 6. We also show the epoch folded history result of the category C(B) source 4U 1907-1-097 
detected in the BATSE EOT light curves by Laycock et al. (2003). (b) Flux histories in the indicated energy 
band for the two category B sources LMC X-4 and SMC X-1 epoch folded at the appropriate superorbital 
periods as given in Table 6. Note the effect of the superorbital period shift in SMC X-1, which smears the 
profile over long times. 

Fig. 16.— Images of selected category B sources: (a) 3C 273 (20-50 keV), an AGN, (b) 4U 1812-12 (Ser 
X-2) (20-50 keV), a LMXB, (c) PSR 1509-58 (50-100 keV), a SNR, and (d) EXO 0748-676 (20-50 keV), a 
LMXB. 

Fig. 17.— Broadband spectra for (a) the LMXB Cir X-1 and (b) the SNR PSR 1509-58 from the JPL EBOP 

(open squares) and MSEC EOT measurements (filled circles) for Tnmcated .Julian Dates 8393-8800. For 
comparison, we also show the 9.1y average 4- band fluxes from Table 4, with systematic errors applied, and 
converted to photon fluxes using the multipliers in Table B2. 
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